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The foundations of the chaotic scattering theory for transport and reaction-rate coefficients for 
classical many-body systems are considered here in some detail. The thermodynamic formalism of 
Sinai, Bowen, and Ruelle is employed to obtain an expression for the escape rate for a phase space 
trajectory of a system to leave a finite open region of phase space for the first time. This expression 
relates the escape rate to the difference between the sum of the positive Lyapunov exponents and 
the Kolmogorov- Sinai entropy for the fractal set of phase space trajectories which are trapped 
forever in the open region. This relation is well known for systems of a few degrees of freedom, and 
is here extended to systems with many degrees of freedom. The formalism is applied to smooth 
hyperbolic systems, to cellular-automata lattice gases, and to hard-sphere systems. In the latter 
case, the geometric constructions of Sinai et al. for billiard systems are used to describe the relevant 
chaotic scattering phenomena. Some applications of this formalism to nonhyperbolic systems are 
also discussed. 
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In a previous paper (I) [1], we extended the chaotic scattering, or escape-rate method of Gaspard and Nicolis [2] for 
"^"the coefficient of diffusion of the moving particle in Lorentz gas systems so as to apply to a wider class of transport 
^^coefficients for a simple, classical fluid and to chemical reaction rates. The line of the argument in I was as follows: (1) 
^-rOne can associate to every transport process in a fluid a microscopic, dynamical quantity called a Helfand moment 
^~^[3], (2) For large enough systems and for long enough times, the mean square fluctuations of the Helfand moments 
ly-^about their initial values grow linearly with the time t, in the case that normal hydrodynamic processes take place in 
^\the fluid. The coefficient of the linear term in t for each Helfand moment is, apart from numerical factors, the relevant 
""^transport coefficient. (3) This "diffusion" of the Helfand moment can be regarded as the result of chaotic scattering 
^rjrocesses in an appropriate phase space. The dynamics of the Helfand moments become more diffusion-like as the 
'^system size increases and as the time becomes longer, due to the occurrence of more and more individual scattering 
I events taking place as the time gets longer. (4) One can characterize this diffusion-like process in phase space in terms 
Q f an escape rate of trajectories from regions where the Helfand moment lies within some prescribed bounds into 
(—regions where the value lies outside these bounds. This is the phase space analog of characterizing Brownian motion 
Qby the rate at which Brownian particles pass into an absorbing boundary. (5) One then relates this escape rate for a 
Helfand moment to a transport coefficient on one hand, and on the other hand, to the sum of the positive Lyapunov 
exponents and to the Kolmogorov-Sinai (KS) entropy that characterizes the set of trajectories in phase space where 
the Helfand moment lies forever within prescribed bounds. These trajectories form an unstable fractal set in phase 
space, a fractal repeller, denoted by , where a denotes the transport coefficient associated with a time dependent 
Helfand moment G[°^ , which on the fractal repeller remains within the bounds 

-X/2 < G[ a) < + X /2 (1) 

where x ls re& l an d positive. 
The main result of I is that 

a = lim (A hm f V Xi(TZ^) - h KS (TZ^)] (2) 

x ' A;>0 

where • >o ^* C^x^ ) ^ s the sum over a ^ positive Lyapunov exponents for the trajectories on the fractal repeller, 
and /iks(7?-x^) ^ s the Kolmogorov-Sinai entropy for trajectories on the repeller. Here the limit V — > oo denotes the 
thermodynamic limit, taken before the limit \ ~~ * 00 • 
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The purpose of this paper is to discuss the chaotic scattering theory which leads to the relation between the escape 
rate 7^™' of the Helfand moment from the region (1), and the dynamical quantities \i(1Z^) and /iks (T^i™' ) , that is 



This result is already well known, due to the work of Kantz and Grassberger [4], Eckmann and Ruelle [5], Bohr and 
Rand [6], Tel and co-workers [7], Grebogi, Ott, and Yorke [8], and Kadanoff and Tang [9], for systems with a few 
degrees of freedom. Our goal here is to extend these previous discussions to systems with many degrees of freedom, and 
to provide the necessary foundations for the application of this method to a number of systems of physical interest. 
These systems include smooth, hyperbolic as well as some nonhyperbolic systems, cellular-automata lattice gases, 
and systems of hard-sphere particles. In the latter case a more delicate analysis of phase space trajectories than that 
needed for smooth systems must be carried out due to the discontinuous nature of the hard-sphere potential. 

Our work will, in the main, draw upon two mathematical developments, both stimulated by problems of interest 
to statistical mechanics. The analysis of the connection between chaotic scattering and escape rates presented here, 
like that of some other workers [6,7], is based upon the thermodynamic formalism for describing dynamical systems 
[10-13]. In this analysis, the escape-rate formula, Eq. (3), is a consequence of the properties of the Ruelle, or 
topological, pressure that arises naturally in the thermodynamic formalism. Our principal contribution here is to 
present the method in a way that it applies to systems with many degrees of freedom, to emphasize some of the 
mathematical literature that has a direct bearing on problems of physical interest, and to show how the method can 
be applied to some systems of current physical interest, including some nonhyperbolic systems. 

The other mathematical development which we will use is the analysis of the dynamical properties of hard-sphere 
systems based on a "ray optics" description given by Sinai, Bunimovich, Chernov and others [14,15]. This description 
involves a careful study of the differential geometry of the phase space trajectories for hard-sphere systems and leads 
to a continued fraction expansion of the so-called second fundamental operator. This result is of deep fundamental 
interest for describing the ergodic properties of hard-sphere systems. The second fundamental operator can be used 
to compute quantities such as Lyapunov exponents and Kolmogorov-Sinai entropies for such systems. An important 
application of these and related methods to the triangular Lorentz gas has recently been given by Gaspard and 
Baras [16,17]. Furthermore, van Beijeren and Dorfman [18] have shown that there is a clear connection between the 
differential geometry analysis of Sinai et al. and the kinetic theory of gases. This leads us to believe that it should be 
possible to reformulate kinetic theory in such a way that its reliance on the chaotic behavior of gases is clearly well 
founded upon a mathematical analysis of the underlying dynamics rather than assumed, as is usually done. This task 
awaits further study and will not be considered here. 

The plan of this paper is as follows. In section II, we present the thermodynamic formalism for smooth hyper- 
bolic systems, derive the escape-rate formula from the properties of the Ruelle pressure, and extend the analysis to 
nonhyperbolic systems. In section III, we apply the geometric methods of Sinai et al. to describe the dynamical 
properties of hard-sphere systems in terms of the second fundamental form. In section IV, we show that the ideas of 
the thermodynamic formalism can easily be applied to cellular-automata lattice gases, which are of interest as simple 
models that exhibit hydrodynamic phenomena typical of real fluids. This topic is explored in considerable detail in 
related work by Ernst, Dorfman, Nix, and Jacobs [19]. In section V we present a brief discussion of related methods 
for describing transport phenomena in terms of the underlying dynamical properties of the system, and we conclude 
in section VI with a number of remarks outlining open problems and directions for further work. 



II. LARGE-DEVIATION OR THERMODYNAMIC FORMALISM FOR DYNAMIC INSTABILITIES 

In this section we apply the large-deviation, or thermodynamic, formalism of Bowen, Ruelle, and Sinai [10-13] to 
study the statistical properties of the phase space trajectories of systems with smooth potential energies and which 
obey Hamiltonian mechanics. The purpose of this formalism is to use statistical properties of the trajectories in order 
to construct invariant probability measures. These measures in turn can be used to compute phase space averages 
which are needed for computing effects of the chaotic scattering phenomena of physical interest here. The large- 
deviation formalism which we use goes beyond a linear or quadratic deviation of some quantity from its reference 
value. For example, the chaotic behavior of a system depends upon exponential separation of trajectories in phase 
space. The large-deviation formalism is designed to treat such circumstances. 




(3) 



A;>0 
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A. Linear Stability and Lyapunov Exponents 



1. The phase space and the tangent space 

Consider a mechanical system of TV particles with Hamiltonian H(q, p) where (q, p) is a 27V/ dimensional vector 
space with q = (q 1; q2, • ••, qjv) where q, is the (/ -dimensional) coordinate of particle i, p = (p 1; p 2 , • • •, pjy) an d Pi 
is the (/-dimensional) momentum of particle %. The Hamiltonian function has the form 

N 2 

^p)=E|+ y w ( 4 ) 

where V(q) is the potential energy of interaction between the particles. To avoid bound states, and orbiting collisions, 
we assume that the potential energy is a sum of central, short range, repulsive pair potentials. Hamilton's equations 
of motion are 

. ag(q,p) . flH(q,p) ... 
q= ^p— ' P = 5q— (5) 

We will suppose either that the boundaries of the system are hard walls with infinite mass whose shape will be 
discussed later, or that periodic boundary conditions are applied. In both cases the total energy is conserved, and in 
the case of periodic boundary conditions, the total momentum is also conserved. We will always consider a system 
at one fixed energy E, and therefore the 27V/ dimensional phase space reduces to the 27V/ — 1 dimensional constant 
energy surface M defined by the condition that H = E. A point on the surface will be denoted by X . The trajectory 
of a phase-space point , initially at Xo, will be indicated by a flow $* such that the phase-space point at a time t later 
is given as Xj=$*Xo. Finally we denote by TAi(X) a linear vector space which is tangent to M at the phase-space 
point X on the constant energy surface. 



2. The fundamental matrix and its decomposition 

An important characterization of the trajectories is given by their stability and, in particular, by their linear stability 
which controls the way infinitesimal perturbations evolve with time in the tangent space. These perturbations can be 
calculated by integration of a coupled set of equations, one for the trajectory that passes through the point X, 

X=F(X) (6) 

which is a simple rewriting of Hamilton's equation, Eq. (5), and one for the a trajectory that deviates by an 
infinitesimal amount from the reference trajectory, X( = <J?'Xo, 

SX = F (X + SX) - F (X) = ^ffi • SX (7) 



to linear order in SX. Since Eq. (7) is linear, all of its solutions can be expressed as 

$*(X 
~dX~ 



SX = $' (X +SX ) - (Xo) = ^ °> ■ SX = M (t, Xo) • SX (8) 



with 



f* dF 

M(t,X ) =T exp^ ^($ r X ) dr (9) 



where a time-ordered exponential is indicated. The fundamental matrix thus obeys the evolution equation 

M(i,Xo) = ^($ i X )-M(i,X ) (10) 

In Eq. (8) we have set X and <5X to be the value of X and SX at t = 0, and M (t,X ) is the fundamental matrix 
solution of Eq. (7). We mention for use in subsection II. D that the construction of the fundamental matrix requires us 
to follow two nearby trajectories through their full evolution over a time interval of duration t. Our primary interest 
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resides in characterizing the rate of separation or of approach of two nearby trajectories in terms of characteristic 
exponents, called Lyapunov exponents. In a stability analysis, stable and unstable directions are determined by the 
signs of the Lyapunov exponents in an obvious way. 

The Multiplicative Ergodic Theorem of Oseledets [20-22] allows us to clearly identify the locally stable and unstable 
directions. This is accomplished through the use of the property of the matrix M(i, Xo) as a multiplicative cocycle. 
That is, M(£, Xo) satisfies the important group relation 

M(t + s,X) = M(i,$ s X)-M(s,X) , (11) 

for any positive time s. This is just a restatement of the time evolution property of the solution of Eq. (7). Next we 
define a Lyapunov homology as a local, linear transformation between the cocycle M and another cocycle m of the 
form 

M(t,X) = C($*X)-m(i,X)-C _1 (X) (12) 
such that the transformation matrices C have no exponential time dependence, that is, such that 

lim -ln||C($'X)|| =0. (13) 

t— >oo t 

One can easily check that the cocycle M satisfies Eq. (11) provided that the cocycle m does also. The purpose of the 
Lyapunov homology is to find conditions under which the cocycle M might be reduced to a diagonal cocycle m such 
that 

M(t,X) = £e fc ($'X)A fc (t,X)f£(X) (14) 

k 

where e^(X) are vectors with components e^^X) = C,&(X) and f& are vectors with components /^(X) = [C _1 (X)]j^. 
Here the superscript T denotes a transpose. These vectors form a set of bi-orthogonal pairs satisfying 

^2 e k = 1 , and fj • e/ = S k i (15) 

k 

The functions Aj,(i, X) are called stretching factors in the directions {e^}. They must also satisfy a multiplicative 
cocycle relation 

A k (t + s,X) = A k (t,$ s X)A k (s,X). (16) 

We emphasize that this construction is not a diagonalization of the matrix M in the usual sense because the vector 
f k is evaluated at the initial point X while the vector e k is evaluated at the final point $*X. It is also important to 
note that e^($'X) and f k (X) are not, in general, mutually orthogonal. The decomposition (14) of the fundamental 
matrix is also referred to as the Mather spectrum [23]. 

When the matrix M is applied to one of the direction vectors one obtains 

M(t,X)-e*(X) = A fc (t,X)e fc ($'X) (17) 

According to the Lyapunov condition Eq. (13), the vector e^($ i X) has no exponential time dependence, so that the 
entire exponential behavior, if any, on the right hand side of Eq. (17) is contained in the function A^ (t, X). Since the 
Lyapunov exponents measure the rate of exponential separation or of approach of two nearby trajectories in different 
directions in phase space, we see that the Lyapunov exponent associated with the direction e k (X) is 

X k (X) = lim iln|Afc(i,X)| (18) 

t—¥00 t 

If the initial perturbation, <5X = e, points toward an arbitrary direction in the tangent space, its time evolution will 
be determined by the largest among the stretching factors for which the scalar products fj-e do not nonvanish as we 
can conclude from the relation (8) defining the fundamental matrix and its decomposition (14). 

In this way Oseledets' theorem shows that the Lyapunov exponent associated with an arbitrary tangent vector e, 

A (X, e) = lim - In ||M (t, X) • e|| (19) 

t— >oo t 
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takes its value from a discrete set called the spectrum of Lyapunov exponents which satisfy 

A« (X) > \W (X) > • • • > A M (X) (20) 

with multiplicities 

mW(X),m< 2 )(X),..,m( r >(X) (21) 
which sum up to the dimension of the tangent space M = dimTA^ = dim.M 

r 

M = J2 m(n) ( x ) = 2N f~ 1 ( 22 ) 

n=l 

The value A^™) (X) is obtained from Eq. (19) when the tangent vector e belongs to the subspace V^"'(X) \ V (™ +1 )(X) 
(that is, the set of points in but not in V'™ +1 '), where {V^™^(X)}^ =1 are nested linear subspaces of the tangent 
space 

T.M(X) = V (1) (X) D V (2) (X) D ■ ■ ■ D V (r) (X) , 

such that m(")(X) = dim[V("H x ) \ V(" +1 )(X)]. The linear subspace V(")(X) is spanned by the set of unit vectors 
{e,(X)} iGJ( „) such that the corresponding Lyapunov exponents in Eq. (20) satisfy A,(X) < A(")(X) for i € I^ n \ 



3. Stable and unstable manifolds 

The stable and unstable manifolds W s (X-o) and W^Xo), respectively, associated with the trajectory at Xo play a 
central role in dynamical systems theory [5,23-25]. These manifolds are defined by 

W s (X ) = {X G M : || $*X - $*X 1| -> for t -> +oc} (23) 

and 

W u (X ) = {X £ M : ||$'X - $'X || -> for t^-oo} (24) 

where ||-|| denotes a distance which we may take to be a Riemannian metric distance on the constant energy surface. 
The stable and unstable manifolds are global objects extending in phase space, as illustrated in Fig. 1. When 
restricted to the vicinity of the trajectory, these manifolds are referred to as the local stable and local unstable 
manifolds. The local stable and unstable manifolds are tangent respectively to the stable and unstable directions 
given by the corresponding vector fields {e^(X)} of the tangent space. The union of all the stable or unstable 
manifolds of all the points of a trajectory U_ 00<t< -|_ 00 W Sj „($'X ) are invariant under the time evolution. 

The concept of stable and unstable manifolds allow us to obtain an alternative method of calculation of the Lyapunov 
exponents. If, in the definition of the global invariant manifold, Eqs. (23) and (24) we replace X by Xo + 5Xo we 
obtain the definitions of the local stable and unstable manifolds in terms of the distance 

llJXtH 2 = ■ M T (i, Xo) • M(i, Xo) • JX 0) (25) 

which defines a positive definite quadratic form. We remark that this quadratic form is strictly internal to the tangent 
space at the only initial point Xo, so that there is no explicit reference to the infinitesimal vector S~Kt at another 
point on the trajectory. As a consequence, the structure of the stable and unstable manifolds, as well as the problem 
of linear stability, can be solved locally at each point X of the flow. Of course, the stability problem requires an 
integration of the linearized trajectory in the tangent plane to the constant energy surface at Xo and this information 
is contained in M T (i, X ) and M(i,X ) in Eq. (25). 

The quadratic form in Eq. (25), M T • M, can be diagonalized in terms of eigenvalues {crj (t, Xo)}, and orthonormal 
eigenvectors {uj (t, X )} as 

2JV/-1 

M T (i,X )-M(i,X ) = u7(i,Xo)<7 i (i,X )u i (t,X ) (26) 

i=l 
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where we have used the fact that the eigenvectors must span the 27V/ — 1 dimensional space TM. (Xo) tangent to 
the constant energy surface at X . The local Lyapunov exponents of the trajectory at initial point X may now be 
defined as 

Ai(X ) = lim IhnTifrXo). (27) 

t— foo Zt 

which gives results which are equivalent to Eq. (18). Depending on the sign of the Lyapunov exponent A, (Xo), the 
corresponding directions are stable (A, < 0), or unstable (A, > 0), or neutral (A, = 0). The neutral directions include 
the direction of the flow (and, if we wished to include it, the direction perpendicular to the energy surface). We 
note that the orthonormal eigenvectors u,(i,Xo) are not directly related to the directions {ejt(Xo)} which are not 
orthogonal. As a result of these considerations, we are able to express the tangent space at X as a direct sum of 
three linear subspaces 

TM (X ) = S u (X ) © S (X ) © £ t (X ) , (28) 

spanned by the unstable, the neutral, and the stable directions, respectively, by merging together the linear subspaces 
V (n) (X ) \ V (n+1) (X ) with A(™> > 0, A(™> = 0, and \^ < 0, respectively. 

4- Local stretching rates 

We notice that the local positive Lyapunov exponents are quantities which are constant on half-trajectories $'X 
with t £ [0, +oo) or (— oo,0] in the sense that A^(X) = Ajt($ s X) for s € M ± , which follows from the cocycle property 
(16) and the definition (18). 

Our purpose is here to introduce the local stretching rates which are underlying the local Lyapunov exponents 
and which could be varying functions along each trajectory. With this aim, we differentiate Eq. (17) with respect 
to time to get an equation of evolution for the stretching factors A^(i,X). Using the evolution equation (10) of the 
fundamental matrix and the bi-orthogonality (15), we obtain 

A*(t,X) =xt(* t X)A /fe (t,X) , or A k (t, X) = exp / X *($ r X)dr (29) 

Jo 

where we introduced the local stretching rate in the direction (X) by 1 

X *(X) = fT(X) — (X)-e,(X) - fT (X ).^(X).F(X) (30) 
These quantities are defined locally at each phase-space point. Accordingly, the local Lyapunov exponents are given 

by 

A*(X) = lim 1 / X *(*'X)dt= lim _L f Xk (&X)dt (31) 

T->oo 1 J T-*oc 21 J_ T 

where we used the time reversibility. We remark that, for linear vector fields, F(X), the first term of Eq. (30) directly 
provides the local Lyapunov exponent A&. The linear stability analysis shows therefore that all of the Lyapunov 
exponents are defined in terms of the local stretching rates Xfc(X). The sum of the local stretching rates 

«(X) = **( X ) ( 32 ) 

\k>0 

defines at each phase-space points X the quantity we call the local dispersion rate. This quantity plays a central 
role in the following considerations. Let us note that some of the local stretching factors may be locally negative 
although the corresponding local Lyapunov exponent is positive so that the sum should extend over all the positive 
local Lyapunov exponents. 

Furthermore, we have the approximate relation 

II V°i ( T > X ) - II A*(T,X) = exp 

<n>l A k >l 

which will be used below. 



/ (** x ) dt 

Jo X k >0 



as T ^ oo. (33) 



1 These local stretching rates do not seem to have already been introduced in the literature. 



6 



B. Symplectic dynamics, the pairing rule of the Lyapunov exponents, and hyperbolicity 

Since Hamiltonian systems are symplectic, the fundamental matrix M(i, X ) obeys the relation 

M T SM = £ 

with the symplectic form £ given by 

* = (-!*)■ <«> 

As a result of this relation, if a is an eigenvalue of M T M then so is a~ l . Accordingly, to each stable direction, there 
corresponds an unstable direction, and vice versa. Their Lyapunov exponents are respectively —A, and +A,. This 
pairing rule also implies that the sum of all of the Lyapunov exponents must be zero for conservative Hamiltonian 
systems and that phase-space volumes are preserved by the dynamics. 

In the case that the Hamiltonian system has periodic orbits, the linear stability of these orbits can be characterized 
by the eigenvalues A, of the fundamental matrix M (T) where T is the period of the orbit. Indeed, the A, are related 
to the preceding eigenvalues by er, = |A,| 2 , so that 

A, = ^lnlA.I . (35) 

It is worth noting that to determine the Lyapunov exponents for non-periodic points Xo, one should, in general 
determine the Lyapunov exponents from the eigenvalues of the symmetric, positive definite form M T (t, X ) M (t, X ) 
but for a point X on a periodic orbit, one may determine the Lyapunov exponents by computing the eigenvalues of 
M (T, X) directly. 

A dynamical system is said to be hyperbolic if all the periodic orbits are unstable of saddle type with non- vanishing 
Lyapunov exponents (except for the zero exponents associated with the direction of flow as well as with any known, 
globally conserved quantities). 2 We shall say that the system is nonhyperbolic if there is a set of periodic orbits for 
which all of the Lyapunov exponents are zero. We will consider systems of both hyperbolic and nonhyperbolic types, 
although there are only a few general statements that can be made about nonhyperbolic systems. 

C. The Large-Deviation Formalism 

The central quantity of interest in the large-deviation, or thermodynamic, formalism is an invariant probability 
measure on the set of phase-space trajectories. Here we outline the construction of the probability measure and define 
a fundamental quantity, the topological pressure. It will be important to keep in mind both the methods and results 
of equilibrium statistical mechanics of classical systems, since this subject motivates a large part of the discussion to 
follow, as was first demonstrated by Sinai, Ruelle, and Bowen [10-13]. In this context, the thermodynamic formalism 
is considered in time instead of space in contrast with equilibrium statistical mechanics. The topological pressure 
plays a role in dynamical systems theory very similar to that of the free energy for statistical mechanical systems. 
In addition it also is essential for establishing an important connection between the invariant probability measure on 
trajectories and the microcanonical measure on the constant energy surface. This connection lies at the heart of the 
large-deviation formalism, which, following the methods and ideas of Sinai, Ruelle, and Bowen, we hope to make clear 
below. 



1. Separated subsets and topological entropy 

We begin by considering the neighborhood of a point X , and suppose that some of the Lyapunov exponents, A, (X ) 
are positive. Then a typical neighborhood of this point of radius e say, will in the course of time be exponentially 
stretched in some directions and exponentially squeezed in others, while the phase-space measure of this small region 



2 We do not include in the definition of hyperbolicity the condition that the stable and unstable directions vary continuously 
with the phase-space point. This condition does not hold for dispersing billiards. 
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remains constant in time. We will want to examine a very small subregion of this neighborhood such that all 
trajectories starting from points in the small subregion will remain within a distance e of the trajectory from X at 
any time over a time interval (— T, +T). Since trajectories separate exponentially rapidly from each other, the size of 
this subregion must be exponentially small, with dimension of order 



e exp 



/ 

J-T 



+T 

t 



Xi(^Xo)^ 



(36) 



where xi is the maximum of the local stretching rates (30). 

We make this argument more precise and generalize it somewhat as follows. We use the notion of (e, T)-separated 
subsets of Bowen and Walters [11,26]. Consider a time T and define a new distance between two points Xi and X 2 , 
Pt(Xi, X 2 ) , on the constant energy surface by 

p T (Xi,X 2 )= max ||$*Xi - $*X 2 || . (37) 

-T<t<+T II H 

If it were to happen that pr(X-i , X 2 ) < e then the trajectories of two phase-space points starting at Xi and X 2 would 
remain within a distance of e over the time interval —T < t < +T. 

We suppose that there exists an invariant set 1Z of interest in the constant energy surface, i. e. a set such that 

^ t (TZ)=TZ (38) 

Consider now a time interval (— T, +T) and a subset of the invariant set, S C 1Z, composed of points which are 
separated by a pr distance larger than e. That is, we construct a set of points S = {Yi, . . . , Y,s} of the invariant set 
1Z such that pt(Yj, Yj) > e for i, j = 1, . . . , S. If the invariant set is compact, i.e. bounded, then one can always find 
a subset S with a finite number of points. This set is called a (e, T)- separated subset of the invariant set. 

Since (e, T)-separated subsets exist with a finite number of points, there exists at least one set, St (e), with the 
maximum number of points, say st(s,TI) points. Since the trajectories separate exponentially with T we expect 
that st (e) will increase exponentially with T also because the points may become separated by smaller and smaller 
distances according to Eq. (36). This rate of growth of St (e) is characterized by the topological entropy defined as 

htov (*) = lim lim sup ^- In s T (e) . (39) 



2. Topological pressure and the dynamical invariant measures 

The idea of the thermodynamic formalism is to introduce a functional of physical observables which is the generating 
functional of the average and of the multitime correlation functions of the given observable ^4(X). The functional is 
called the Ruelle topological pressure and is defined as 

V(A) = lim lim sup — j- In Zt (A, e) . (40) 

with the partition functional 

r+T 



A(^Y)dt 



(41) 



where S is a (e, T)-separated subset of the invariant set 1Z. The topological pressure has remarkable properties. In 
particular, it is a convex functional of the observable, i. e. V\yA + (1 — v)B] > vT(A) + (1 — v)V(B) for < v < 1 
and any two observables A and B. 

When the observable is everywhere vanishing, A = 0, the topological pressure reduces to the topological entropy 
because 

s T (e) = Z T (A = 0,e) , so that h top ($) = V(0) (42) 
Let us now consider another observable B of the dynamical system. The average of this observable is defined by 
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Ha{B) = —V(A + vB) 



v=0 



(43) 



Introducing the definition of the pressure, we get that ha(B) = J B(X)[iA(dX) with the measure 

exp f +T A(&Y)dt 1 f+ T 
» A (dX) = lim lim sup Sup 5 £ J ~ T ' — S(X- &Y)dt dX 

e^o T ^oo ^ s Z T (A,e) 2T J_ T 



(44) 



which is referred to as a dynamical measure. We observe that each trajectory of the (e, T)-separated subset is weighted 
by a Boltzmann-type probability given by 



exp/+jA($ i Y)di 

Zr(A,e) 



7r A (Y,T,e) 



(45) 



where A(&X)dt plays the role of —/3E. Accordingly, such measures have been called Gibbs canonical measure. 
Here, the role of the energy is played by the average of the observable A over the time interval (— T, +T). We 
emphasize that the Gibbs measure, Eq. (44) is determined by a time interval (— T, +T) rather than by a number of 
particles, interaction range, etc., as in equilibrium statistical mechanics. In the limit T^oo the probability measures, 
Eq. (44), tend to measures which are invariant under time evolution, known as "equilibrium states" as proved by 
Bowen and Ruelle for Axiom- A hyperbolic systems [12]. In this limit a connection can be established between the 
measures Eq. (44) and the invariant probability measure on the invariant set, a connection which will be important 
for us subsequently. 

In this formalism, correlation functions between two observables B\ and B 2 can be obtained as second derivatives 



a 2 



■V{A + v 1 B 1 +v 2 B 2 ) 



/+oo 
-oo 



Va(BiB 2 o $') - i_ia(Bi)[ia(B 2 ) 



dt 



(46) 



in a straightforward notation, provided the various limits exist. Moreover, the Kolmogorov-Sinai entropy per unit 
time of the dynamical system with respect to the invariant measure [ia is defined by [24,27] 



hush* a) = limlimsup ( —^=) Sup 5 Tr A (Y,T,e) ln7r^(Y,T,e) 



(47) 



We can now combine Eqs. (43) and (47) with B = A to obtain the fundamental identity of importance to us, 

hKs{HA) = -HA{A)+V(A) (48) 



3. Pressure functions based on the Lyapunov exponents 

By varying the observable A, we can obtain many invariant measures, and it is not yet clear which among these 
many possible measures is the one appropriate to a specific numerical simulation of the dynamical system starting 
from a given statistical ensemble. The answer to this question will be delayed till the next two sections where the 
question is answered for closed and open hyperbolic systems. Nevertheless, one thing that is common to the closed 
and open cases is the special role played by the sum over the local stretching rates multiplied by a parameter — f3 in 
order to emphasize the formal analogy with the Gibbs states 

A(X) = -pu(X) = ~PJ2 »( X ) ( 49 ) 

A;>0 

in terms of the local dispersion rate (32). This observable measures the dispersion of trajectories emanating from 
points in this region over a time interval (— T, +T). The larger the dynamical instability on the trajectory, the smaller 
the probability to visit the neighborhood of this trajectory. This reasoning is at the basis of the choice (49). In this 
case, the pressure functional becomes the following pressure function 

p(/3) = r[-/3Y / x i m} (50) 

A;>0 
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which defines an invariant probability measure up depending on the parameter /3. The sum of averaged stretching 
rates is then given by 

>=»(E«') < 51 > 

A;>0 

Using the time invariance of the dynamical measure fip applied to Eq. (31), we obtain that the averages of the local 
stretching rates are identical with the averages of the corresponding local positive Lyapunov exponent 

VpiXi) = Vp(^i) (52) 
Therefore, the fundamental identity (48) becomes here 

h K s (pp) = P ^ ( A + P W = P M«) + P(P) (53) 

Ai>0 

We notice that the parameter fi has nothing directly to do with the inverse of a thermodynamic temperature. It is 
worth noting, however, that the thermodynamic formalism does make possible deep connections between the calcu- 
lation of thermodynamic functions and dynamical properties, a subject which is discussed at some length elsewhere 
[10-13,28-30]. 

The individual average Lyapunov exponents can be obtained by defining a multivariate pressure function 
P(Pi, • ■ • , Pl) from the observable 

A(X) = - £ ft Xi (X) (54) 

A;>0 

depending on the L = Nf — 1 parameters /3 = (/3i, . . . , The averages Lyapunov exponents are now given by 

-|^(/3) =/i/3 (A,) (55) 

When all the parameters are equal f3 = f3i = . . . = 0l, the pressure function (50) is recovered: P(f$) = P(Pi, ■ ■ ■ 
(see Fig. 2). 



4- Entropy function and Legendre transform 

It is convenient to introduce also a entropy function by considering the number of points of the (e, T)-separated 
subset such that the time average of their associated local stretching rates, calculated over a time interval (— T, +T), 
take values in the interval ((fi, Lpi + d(fi) according to 

NumberjYGcS: ^ Xi(**Y)di G ebpt), z = 1,...,lJ = p(T,S,cp) exp[2T5(<p)] d L <p (56) 

in the limit where e — > and T — > oo, where p(T,S,(p) is a slowly varying function of the time T. The entropy 
function S((p) is known to be a concave function [6,31,32]. 

The entropy function is related to the pressure function P{[3) by a Legendre transform. Indeed, the sum over all 
the points of the (e, T)-separated subset S can be replaced by an integral over (p 

P(f3) = hm -L In / d L tp p(T, S, <p) exp[2TS(cp)} exp(-2T ¥ > • 0) (57) 

T->oo 21 J 

In the limit T^oo, the integral can be evaluated by the steepest-descent method which selects the maximum (pp of 
the function in the argument of the exponential as solution of 



dtp 

Therefore, the pressure function is given by 
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P((3) = S(<pp)-cpp-(3 (59) 
Reciprocally, once the pressure function is known, the entropy function is obtained as 



dP 

S(<p)=P(0 v ,) + <p-0 v with v = - (fi^) (60) 



In particular, the topological entropy is given by 



dP dS 

hto P (*) = P(0) = S(<p ) with Vo = __(o) or ^(Vo)=0 (61) 
and the KS entropy by 

husivp) = S{ipp) since ipp=fip{\) (62) 

which justifies the name. Here also, we can define a univariate entropy function S(<p) = S(<p, . . . , <p), which is related 
to the univariate pressure function by a Legendre transform. 

We need a proper interpretation of the different terms appearing in the fundamental identity (53) , especially of the 
pressure P (f3) and we need to fix the value of the parameter (3. 



D. Closed Hyperbolic Systems 

1. The microcanonical measure as a Sinai-Bowen-Ruelle dynamical measure 

In a closed, conservative, ergodic system the natural measure is the microcanonical measure [i. The appropriate 
invariant measure on this surface is, of course, the microcanonical measure /i eq (cZX) given by 

^<"> = ii=§i < 63 » 

where dcr(X) is a surface area element about the point X and H is the Hamiltonian of the system. The gradient is 
a 2Nf dimensional gradient evaluated on the surface H = E. Our purpose is to identify the microcanonical measure 
with one of the dynamical measures introduced in subsection II. B. With this purpose, we consider one of the small 
regions surrounding a point Y in a finite (e, T)-separated subset S and the probability of this region can be computed 
using Eq. (63) once the region is identified. In subsection II. B, we considered that this same region is a domain on 
the constant energy surface in which the trajectories of all points will remain separated by a distance less than e over 
a time interval (— T, +T), which is known as a ball Bt(Y,s). The microcanonical probability of such a ball can be 
estimated using the results of subsection II. A as follows 

//eq[»T(Y,e)] = /i eq {x £ M : ||**X - $'Y|| < e, V* £ [-T,+T]} 

~ ^qjxeX: ||M(t,Y) -(X-Y) || <e, Vie [-T,+T]} 
= //, q {xeX: £>i(*,Y)|in(t,Y).(X-Y)| 2 <e 2 , Vte[-T,+T]} 



; (+T,Y)>l v / (7;(+Ti y) Ha < (-T,Y)>1 ^(.t.Y) 



exp - J^Zx^oXii&Y^t 



~ ex P " J-T u (® tY ) dt 

~ i"«(Y,T,e) (64) 

At the first line, the definition (37) of the distance has been used. We supposed that e is small enough and we used Eq. 
(25) to get the second line. The third line results from the spectral decomposition (26). The fourth line is based on the 
facts that the quadratic form defines a small ellipsoid with axes determined by the quantities er,(i,Y). Half of these 
quantities increase exponentially for t > and the other half increase exponentially for t < 0, while the vectors u,(i, Y) 
are slowly varying functions of time. The fifth line is a consequence of the relation (27) of the eigenvalues Oi(t, Y) to 
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the local stretching rates and of the pairing rule that, to every stretching rate, there corresponds a contracting rate 
with the same absolute value by time reversibility. Finally, the last two lines result from the definition (49) of the 
observable u(X) and of the probability (45). 

To get this last result, we used the fact that the normalization factor Zt(u,e) is a slowly varying (subexponential) 
function of time in the case of closed hyperbolic systems. Indeed, the sum of all the probabilities (64) is approximately 
constant. On the other hand, all the equalities in Eq. (64) hold up to factors which are slowly varying with time. 
According to the fifth line, the sum of the exponential factors involving the local stretching rates is slowly varying. 
With the last line, this implies that the dynamical partition function Zt(u,e) is also slowly varying with time. 
Therefore, all of the exponential dependence should be completely taken into account by the dispersion factor given 
on the right-hand side of Eq. (64). This observation suggests that the member of the family of measures, Eq. (44) 
based on the observable A = —u which corresponds to the value (3 = 1 is the natural invariant measure corresponding 
to the microcanonical measure 



(65) 



This result is general and extends to open hyperbolic systems. 

We now develop some consequences of this identification. Since the system is closed and the equilibrium measure 
of the constant energy surface is finite, both the microcanonical and the dynamic measures must be normalizable to 
unity, say. As a consequence, we already concluded that the normalization factor Zt(u,e) does not depend on time 
in an exponential way 3 , but increases in a subexponential way with T when (3 = 1, and, as a result the pressure at 
this value of (3 must be zero, 



P(l) = for closed systems 
Inserting this result in the fundamental identity, Eq. (53) with (3 = 1 we recover Pesin's identity [5,33] 

/iKs(^eq) = ^2 /W^i) 
A;>0 



(66) 
(67) 



so that the KS entropy of the invariant natural measure is the sum of the positive Lyapunov exponents averaged 
over the same measure. We remark that the dynamical invariant measure which is the microcanonical measure is 
absolutely continuous with respect to the Lebesgue measure along the unstable manifolds. When this property holds 
which is a corollary of the Pesin identity, the dynamical measure is referred to as a Sinai-Bowen-Ruelle (SRB) measure 
[5,30]. 

2. The pressure function for closed systems 

Another remarkable and useful identity for the pressure P(f3) can be obtained in an alternative way as an average 
over the microcanonical measure (63) according to 
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P(/3) = lim sup — In / fi eq (dX.) exp 



T->oo 2T j M 



(1-/3) f +T Xi(**X)dt 

" T Xi>0 



(68) 



The original definition of the pressure can be recovered as follows. Considering an (e, T)-separated subset S of the 
phase space Ai, the integral can be discretized into a sum over the points {Y} of S replacing the volume elements dX 
by small balls Bt(Y,e). This sum would be equal to the integral after the appropriate limits are taken. Therefore, 
the right-hand side of Eq. (68) is given by 



lim lim sup -^ln Sup 5 V ^ eq [B T (Y,e)] exp [(1 - /3) [ + u(^Y)dt 
E-s-o j,^^ 2T ^-? s [ J-T 



(69) 



Now, according to Eq. (64), the probabilities of the small balls are exponentially growing like exp — 

which introduce an extra inverse power of the dispersing factor and explains the power (1 — (3) in Eq. (68). Finally, 

we obtain the definition of the pressure, namely 



3 As it would if particles were escaping from the system at an exponential rate. See subsection II. D. 
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P((3) = lim limsup r 

e^O y^oo 2T 



In Su P< 5 ex P -Z 3 / w(*'Y)di 
Yes L , '~ T 



(70) 



We note that for the purpose of numerical calculations, we can substitute the stretching factors <jj(i, X) for the 
exponentials of the integrated local stretching factors and use 



P(f3) = lim ±ln / /x eq (dX) JJ v^P) 1 



(71) 



instead of Eq. (68). Here, there is no factor of two in the denominator because the stretching factors correspond to 
the time interval (0, +T) rather than (— T, +T) as before. 

Equations (68) or (71) provide a convenient way to compute the pressure function as a phase-space average of an 
expression involving the local stretching factors. A typical pressure function is illustrated in Fig. 3a for a closed 
hyperbolic system. Note that it is a convex, monotonic function of (3 and vanishes at (3 = 1. This is not true of 
open systems. The KS entropy is obtained by finding the slope of the pressure at (3 = 1, and carrying out the linear 
extrapolation illustrated in Fig. 3a. Finally, we note that the value of P(f3) at f3 = 0, is the topological pressure 
^top(^) of the dynamical system, as can be seen from Eqs. (39), and (61). The fact that /iks < h top puts a bound on 
the slope of P{(3) at (3 = 1 



3. Generating functions of transport coefficients 



In paper I, we showed that each transport or rate coefficient a is associated to a diffusive-type motion for a 
corresponding Helfand moment . Within the large-deviation formalism, it is possible to characterize in detail the 
random time evolution of these Helfand moments by introducing, the generating function 



Q(k,/3) = r lim ^ ln^jexp^G^-G^)]} 



(72) 



In particular, if there is no drift of the moment, i. e., if d K Q\ K _ Q = 0, the transport coefficients with respect to the 



SRB measure fip is given by 



18 2 Q 



2 8k 2 



K=0 



(73) 



The standard transport coefficient is given by the value at (3 = 1 at which we recover the microcanonical measure. 
The higher derivatives give the higher-order transport coefficients like the Burnett coefficients [34] . 
We can also define a large-deviation function according to 



f G (a) - G (a) 

W Ye5: +T 2 T ~ T t^ri + dr,) 



p(T,S,r,,P) exp^Ttffa./J)] d V 



(74) 



in the limit T — > oo. In Eq. (74), p is a slowly varying function and the Helfand moment is calculated for the 
trajectory Y of an (e, T)-separated subset of the phase space Ai. The function H(r),f3) is related to the generating 
function by a Legendre transform according to 



8H , 8Q 
— and 77 = — 

8rj 8k 



Q{k,/3) = H(t],0) + k T] , for k = 
These considerations are particularly useful in the study of anomalous transport [35] 



(75) 



E. Open Hyperbolic Systems 

1. Escape-time functions and the repeller 



The main focus of our analysis in Paper I was on the chaotic behavior of open systems, i.e. systems with trajectories 
which escape from a bounded region of phase space in a finite amount of time. The behavior of such systems can be 
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easily visualized by considering an escape-time function as shown in Fig. 4 for a simple model of a logistic map with 
escape from the interval < x < 1. The connection between transport coefficients and dynamical quantities is based 
on the escape-rate formula, Eq. (3) , and here we describe the derivation as a special case of the fundamental identity, 
Eq. (48). 

An open system is defined by a trajectory dynamics taking place in a phase-space region B, with a boundary dB. 
We consider initial conditions X , within B, and define an escape time T B + ^(X ) as the first time when the trajectory 
crosses the boundary 

$*Xo G Bfor < t < T B +) (X ) (76) 
$ T< +) (x ) Xo g gg^ and 
$ t<+>(x ) + * Xo ^ dB,B 

for arbitrarily small S > 0. Since we are considering time-reversible, Hamiltonian systems, we may also define a 
"negative" escape time for an initial condition X € B in an analogous way to that in Eq. (76) by 

$* (X ) G B for T B _) (X ) < t < (77) 

$ t b _) ( x o)(x ) e dB, and 
$ t<->(x )+* (Xo) 

for arbitrarily small S < 0. The escape time is a highly singular function of the initial conditions if there exist periodic 
and non-periodic trajectories which are forever trapped inside the domain B such that $'X € B for — oo < t < oo. 
Since most of the trajectories are expected to escape from B for hyperbolic systems, these trapped set of trajectories 
form a set of zero Lebesgue measure. This set - called the repeller - may contain a subset of chaotic trajectories in 
which case the repeller is a fractal set. Because of hyperbolicity, the points on the trapped trajectories must have 
stable and unstable manifolds. All of the initial conditions in B belonging to the stable manifolds of the trapped 
orbits remain in B for all positive times. For these initial conditions the escape time is infinite. However these initial 
points form a set of zero Lebesgue measure so that the escape-time function takes finite values for almost all of the 
points in B (see Fig. 4). 

We now consider the set of initial conditions for which the escape time is larger than a predetermined time T > 0, 

T ( +\T) = {x eB: T < T B +) (X )} . (78) 
This set contains all of the members of the statistical ensemble which are still inside the domain B at the time T, i.e. 

ri +) (T)= P| $"*(£). (79) 

0<t<T 

In the long time limit, this set contains the trapped trajectories of the fractal repeller 1Zb and their stable manifolds, 
restricted to the domain B: 

lim T ( B + \T) = C\[W s (Tl B )}[]B (80) 

where Cl[-] denotes the closure of the set. As a consequence, the set T^ +) (T) under goes a fragmentation into smaller 
and smaller sets as time increases to end up as the fractal set (80). In a similar way, we can define a set of initial 
conditions such that the time reversed motion has a "negative" escape time which satisfies |T B ^(Xo)| > T 

T { B \T) = {X„6B: IT^Xq)! > T, < o} = p| *-*(B) (81) 

-T<i<0 

In analogy with Eq. (80), we have that 

lim T ( -\T) = C\[W u (n B )]f]B (82) 

The intersection of the two sets, 
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T B (T) =T<" ) (T)nT< +) (T) = p (83) 

-T<i<+T 

contains all of the members of the statistical ensemble which are inside the region B over the time interval — T <t< 
+T. The repeller is defined as the set of points Tg(T) asT^oo 

lim T B {T) = C\[R B ] f)B (84) 

T->oo 1 1 

It is instructive to consider this construction for simple systems such as a two-dimensional Smale horseshoe map of 
the unit square as discussed by Lanford, Tel, and others [25,36,37]. One sees that for this map the sets T B +) (T) 

and Tg \t) are thin strips parallel to the stable and unstable directions, respectively, each becoming the product 
of a Cantor set with a one-dimensional interval as T — > oo. The repeller is the intersection of these two sets, and is 
a Cantor set which can be coded as a bi-infinite sequence of zeros and ones. The dynamics on the repeller is then 
isomorphic to the left-Bernoulli shift on these bi-infinite sequences. 



2. The nonequilibrium invariant measure of the repeller 

We now define the invariant probability measure of the repeller. We first consider a probability measure ^o(dX) on 
the region B corresponding to a particular statistical ensemble of No initial conditions at phase-space points {X^)} 
where these points are distributed uniformly in B with respect to the microcanonical ensemble, say. This measure vq 
is taken to be of the form 

u (dX) = lim — V 6 fx - X (j) ) dX (85) 

JVo-voo Nn \ / 

Of these N initial points, the number Nt still contained in B over the time interval (0, +T) decays according to 

lim ?f = ,„[T< + »(r)] = /' n,(dX). (86) 

We note that the limit N — > oo is essential to define a smooth function of time T, since Nt typically shows large 
statistical fluctuations when Nt < 10, and drops to zero after a finite time. Similar expressions hold for the time 
intervals (-T,0) and (-T,+T). 

Assuming that almost all the trajectories escape, i. e. that lim^oo fo[Tg"' (T)] = 0, the decay curve, Eq. (86), 
may be exponential or slower than exponential in general systems. However, in hyperbolic systems (where all orbits 
are of saddle type) the decay is exponential. Thus we can define an escape rate according to 



7= 



lim ^ In „ [X W (T) n T ^ (T)] (87) 
limjlni* [r B +) (T)] (88) 
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= - T lim o ^lni^[T^ ) (T)j (89) 

All these definitions are equivalent because of time reversibility. For noninvertible systems like ID maps, only Eq. 
(88) is of application. 

In the long-time limit, the trajectories remaining in the domain B are distributed according to a probability measure 
which is invariant for the dynamics on the repeller. In order to construct such an invariant probability measure on the 
repeller, we use ideas familiar from ergodic theory. In the usual arguments of ergodic theory one considers the time 
average of some dynamical quantity for a system whose phase-space point is confined to a constant energy surface. 

If the system is ergodic then the long time average of any dynamical quantity is equal to its ensemble average taken 
with respect to the microcanonical measure. Consider now a system whose phase-space trajectory is confined to the 
repeller. If the trajectories on the repeller are ergodic with respect to a natural measure on the repeller, then the long 
time average of any dynamical quantity on the repeller should be equal to the ensemble average of this quantity with 
respect to the natural measure. To construct the natural nonequilibrium measure on the repeller, then, we begin with 
the definition of the time average of some observable A(X.) on the repeller as 
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N'T j rp 

(lne (A)=\im \un^J2^[ M^ U) )dt. (90) 

T->oo -ZVt-»oo 1\x 21 J _ T 

In Eq. (90), the sum extends over the Nt phase-space points whose trajectories remain in region B over the time 
interval — T <t< +T. The time average can be rewritten in terms of the initial measure as 

Hn e (A) = J A(X)^ e (dX) (91) 

with 

/ine(dX) = lim 1 f MdZ) /t b( t)(Z) ^= f +T S(X - &Z)dt (92) 

T-s-oo i/ [I JJ J 11 J_ T 

Here Ir(X) is the characteristic function of the set T in phase space, i.e., Iy(X.) = 1 if X € T and = 0, otherwise. In 
this way, we have defined a normalized probability measure on the repeller /t ne . This is the desired invariant measure 
with the fractal repeller as its support. We also note that the subscript 'ne' of /x ne refers to the nonequilibrium 
character of the natural invariant measure. This measure is the natural generalization of the microcanonical canonical 
ensemble measure to open systems where the dynamics takes place on the fractal repeller, Using this measure, 
we can consider the long time limit and define the average of an observable according to Eq. (91). This procedure 
amounts to performing the statistics on the set of Nt initial conditions which are still in the domain B at time T. As 
T^oo these trajectories approach more and more closely to the trajectories on the repeller. As a consequence, /i ne 
is an invariant probability measure on the repeller. We will return to this invariant measure presently. 



3. Connection with the dynamical invariant measure and the pressure function 

Now, we wish to construct a dynamical measure on the repeller similar to that used for closed systems in subsection 
II. C. This will enable us to continue the development of the previous sections so as to the fundamental identity, Eq. 
(48) to open systems. The thermodynamic formalism suggests to consider the family of dynamical invariant measures 
defined in analogy with Eq. (44) 

exp —3 [ +T u(&Y)dt i f +T 
^(dX) = lim lim sup Su P<s £ ±f-jr — — S(X - &Y)dt dX (93) 

with the dynamical partition function 

Z T (/3,e) = Su P< 5 ex P -I 3 j u(&Y)dt (94) 
Yes " T 

in terms of the dispersion rate u = 2)2\->oXi- Let us emphasize that the local stretching rates as well as the local 
Lyapunov exponents are well defined for the trajectories of any subset S of the repeller TZb since those trajectories 
remain in the compact domain B forever. A disadvantage of this definition is that the repeller 1Z needs to be known 
since a subset S of it is considered. Physically and numerically, the repeller of systems of scattering type appears out 
of the dynamics after a very long time as above described. The previously defined invariant measure avoids this a 
priori knowledge of the repeller since the invariant measure (92) is automatically constructed by the time evolution 
We now show that both measures are equivalent for open hyperbolic systems when (5 = 1. We do this by using the 
neighborhoods of the points of a (e, T)-separated subset, S, of the repeller 1Z. Note that we take the time T to be the 
same in both sets so that we can closely approximate the repeller by means of the separated subsets as we take the 
limit T — > oo. In this situation, we can transform the expression (92) into 

^(dX) = lim lim sup Sup, £ ^I^A n T B (T)] 1 f +T #l 

It is clear that this measure is normalized. Since the points Y belongs to a (e, T)-separated subset of the repeller TZb, 
they belong to the repeller itself. If the points Y are not close to the boundary dB and if e is small enough, the balls 
Bt(Y,£) are contained inside Tb(T) so that 
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i^[BT(Y,e)nT B (T)]~wb[B T (Y,e)]~exp- / + u(**Y)dt (96) 

i-T 

where considerations similar to the ones of Eq. (64) have been here applied to vq. As a consequence, Eq. (95) becomes 

„ exp- f + ^M($'Y)di l f+ T 

/ine(dX) = hm limsupSups V -^fr-s / - **Y)dt = ^(dX) (97) 

£->o Z T (l,e) 2T J_ T 

Accordingly, we recover the dynamical invariant measure corresponding to the value (3 = 1 as in the case of closed 
hyperbolic systems [see Eq.(65)]. 

We note that the normalization factor is here required for the following reason. In order to estimate the normal- 
ization factor, we decompose the set Tg(T) into small balls Bt(Y,e) centered on the points Y of an (e, T)-separated 
subset S as 

H)[T B (T)] ~ Sup 5 Y, »o[Bt(Y,s) n T B (T)] 

YeS 

r+T 

~ Su P<s ^ exp - / w($*Y)di = Z T (1, e) ~ exp 2TP(1) (98) 
Yes , ' _T 

Comparing with the definition Eq. (87) of the escape rate 7, this result shows that the normalization factor decays 
exponentially like Zt(1,e) ~ exp(— 2T7) and, moreover, that the pressure function P(f3) has, for an open system of 
hyperbolic type, the value 

P(l) = -7. (99) 

Moreover, we can now define the important quantities of the thermodynamic formalism of subsection II. A and the 
fundamental identity (53) follows once again. We can now identify all of the quantities appearing in this equation at 
3 = 1. For this case we find 

7= Y A'ne(Aj) - h KS (ll„e) ■ (100) 

A;>0 

Eq. (100) is the escape-rate formula for an open system, giving the relation between the escape rate, 7, and the 
sum of the positive Lyapunov exponents and Kolmogorov-Sinai entropy for trajectories on the fractal repeller, TZb, 
using the natural measure \i\ = /i ne . We remark that we used the notations \i(TZts) = [i n e(\) = Mi(^i) and 
/iks(^-b) = /iKs(^ne) = ^ks(^i) in Eq. (2). This result generalizes Pesin's formula to open hyperbolic systems. 
Indeed, when the system is closed, the escape rate vanishes 7 = and Pesin's formula is recovered as well as Eq. 
(66). Here, the average of the sum of positive Lyapunov exponents over the nonequilibrium invariant measure can be 
calculated as before as a derivative of the pressure function 

E = ^m 1 (ion 

A;>0 

Let us remark that the non-equilibrium invariant measure is not absolutely continuous with respect to the Lebesgue 
measure along the unstable manifolds but singular with a fractal for support. In this way, the nonequilibrium invariant 
measure differs from the equilibrium one. 

Let us now describe a pratical way to calculate the pressure function in open systems. As we have seen, the pressure 
function require the knowldege of the local Lyapunov exponents. In practice, we only know the stretching factors 
(Ti(T, X( J ') of the trajectories of the initial ensemble (85). In analogy with Eq. (71), we propose here 

P((3) = lim 1 In / JJ v^X) '"^(dX) (102) 

This definition is equivalent to the previous one for the following reasons. By time reversibility, we can convert the 
average over the forward set T B + '(T) into an average over the set Tg(T) by considering the time interval (— T, +T) 
rather than (0,+T). Using an (e, T)-separated subset S to dicretize the integral, the right-hand member of (102) 
becomes 
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limlimsup-^lnSups V v [B T (Y,e) nT B (T)]exp(l -(3) [ u($*Y)di (103) 
s^o j,^^ 2T ^ s ' J_ T 

where we used the estimation (33). According to Eq. (96), we recover the original definition of the pressure 

P((3) = lim lim sup -L In Z T ((3, e) (104) 

with the partition function (94). 

F. Nonhyperbolic Systems 

In many cases the system possesses periodic orbits with vanishing, but non trivial Lyapunov exponents. This is 
the case, in particular, for the stadium billiard, the Sinai billiard, and the hard-sphere gas with periodic boundary 
conditions or placed in a rectangular box. The periodic orbits that we have in mind are those special orbits where 
the particles move without collisions or bounce between parallel walls, for example, as illustrated in Fig. 5. These 
special orbits form sets of zero Lebesgue measure so that they do not prevent the system from being ergodic, or 
a K-flow, or from having positive average Lyapunov exponents. However, all the Lyapunov exponents of these 
special periodic orbits vanish because perturbed trajectories may separate from the periodic orbit in an algebraic 
way: ||$*(X P . .-|-<$X) — $*(X P . .)|| ~ t. As a consequence, if the system is open and the domain B contains such a 
marginally unstable periodic orbit, the decay in Eq. (86) is nonexponential as 

MT ( B + \t)} ~ (105) 

where Nf is the total number of degrees of freedom of the system, and the escape rate vanishes [17,38,39]. (We 
note that this algebraic decay results from purely geometric effects - such as appear when particles travel down long 
corridors without hitting anything- and thus has no direct connection to the long-time tails which appear in the time 
correlation function expressions for transport coefficients [40,41].) 

The large-deviation formalism is very useful in this context because the pressure function is still non-trivial in this 
case. If the Lyapunov exponents vanish in some regions of phase space, some of the dispersion rates u(Y) may be 
vanishing for a few values of Y in the (e, T)-separated subsets. For large and positive values of (3, the partition 
function (94) is dominated by the few terms with u(Y) = because all the other terms are exponentially vanishing 
in the limit T — > oo. As a result, the pressure function is equal to zero for large values of f3. On the other hand, for 
negative values of (3, the terms with nonvanishing dispersion rates u(Y) dominate the sum in the partition function 
so that the pressure is then positive and nontrivial. There exists a critical, lowest value of (3 above which the pressure 
vanishes. 

P(f3) = for f3 c < 13 (106) 

Schematic pressure functions are depicted in Figs. 2c, 2d for nonhyperbolic closed and open systems. For closed 
systems, we observe that the pressure is zero above (3 = 1 which is in agreement with the vanishing of the escape rate 
due to Eq. (105), so that the critical value is f3 c = 1. Moreover, the sum of positive Lyapunov exponents must be 
defined as the left-sided derivative of the pressure at (3 = 1 in the case of a closed system. 

The discontinuity in the shape of the pressure function is the evidence of a phenomenon of dynamical phase 
transitions - so called in analogy with statistical mechanics [7,28-30]. We have the following interpretation. The 
system is described by the continuous family of invariant measures, [ip. The parameter (3 acts like a filtering parameter. 
When (3 > (3 C , the measure is concentrated on the regular trajectories which have vanishing Lyapunov exponents. For 
these measures, we can talk about an ordered or regular phase since the corresponding invariant states are ordered. 
On the other hand, when (3 < (3 C , the measure up gives dominant probability weights to the nonperiodic trajectories 
which are uncountable. In this case, we can speak of a chaotic or a disordered phase. Tuning the parameter (3 therefore 
reveals the chaotic features of the dynamics. This is currently done when we refer to the topological pressure per unit 
time, i.e., P((3 = 0), as an indicator of chaos. 

However, if the pressure function is analytic away from critical points, it has the further advantage that it can be 
extrapolated from below criticality up to the value at (3 = 1 so as to define a supercritical measure, for instance, 
in the case of open nonhyperbolic systems. In this way, it is possible to define an effective escape rate, 7 e ff, as well 
as an effective value for the sum of positive Lyapunov exponents at (3 = 1. The effective rate can be evidenced in 
numerical simulations from the transient behavior of the decay function Eq. (86) . In many nonhyperbolic systems 
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with a large number of degrees of freedom, the power law decay Eq. (105) may remain a very small effect which 
is visible only after extremely long times because marginally unstable periodic orbits are very rare and the decay 
could appear exponential within statistical errors. In nonhyperbolic systems, although we may know from theoretical 
arguments that the escape rate defined by Eq. (87) actually vanishes, the concept of an effective escape rate is useful 
to characterize in a rigorous way a numerically observed exponential decay. 



III. THE HARD-SPHERE GAS 



In this section we will outline the methods that must be used to apply the considerations of section II to a gas of 
hard spheres. The treatment of the hard-sphere gas as a dynamical system requires a different analysis than that given 
in section II for smooth dynamical systems or in section IV for stochastic dynamical systems. Sinai, Bunimovich, 
Chernov and others [14,15] have developed a very beautiful description of billiard systems in order to analyze their 
ergodic properties. Here we will use the methods of Sinai et al. in order accomplish two goals: (1) To present the 
Sinai method in a somewhat elementary way in order to acquaint a more physically oriented audience with these 
techniques which we believe are very useful for the description of the dynamics of hard-sphere systems; and (2) To 
show how the sum of the positive Lyapunov exponents can be expressed in terms of the average value of the trace 
of a certain matrix, called the second fundamental operator or curvature matrix. With this result one can begin to 
apply the escape-rate method to hard-sphere systems. 

In this section we will discuss the theory of billiard systems from the point of view of geometrical optics by studying 
ray trajectories and perturbations thereof. Various notions of elementary differential geometry will be important here, 
and the reader is advised to consult [42] for an introduction to this subject. We also mention a recent, and extensive 
review of billiard systems by Tabachnikov [43] . 



A. Definition of the billiard 



The hard-sphere gas can be thought of as a generalized billiard in a configuration space of dimension Nf where 
N is the number of hard spheres, each of mass m and diameter d, and / is the spatial dimension of each "sphere", 
i.e. / = 2 for disks, / = 3 for spheres, etc. The positions and velocities of all of the spheres are given by the Nf 
dimensional vectors q, and u, 

q= ( r i, r 2, r 3,-, r jv) (107) 

and 

u = ( v l, v 2, v 3>-, v jv) (I 08 ) 

respectively, where r,, and v, are the position and velocity of particle i. The Hamiltonian for the system is 

N 1 

-mv 2 . (109) 

i=\ 



The hard-sphere condition requires that 



>d for j (110) 



and that the particles undergo instantaneous collisions when |r,— r j- 1 = d. Since energy is conserved in elastic collisions, 
the magnitude of the velocity u remains unchanged in a collision. We can then rescale the velocity to a unit value 
u 2 = 1. Then up to a rescaling factor, the dynamics is the same on every energy shell, H = E. 

The positions at which the particles collide |r,— r,-| = d defines a hypersurface dQ in the Nf dimensional space 
which is the border of the billiard Q. We denote by n(q) a unit vector which is perpendicular to the hypersurface dQ 
at the point of impact q and which is directed inside the billiard Q. If er(q) = is the equation of the hypersurface 
dQ, the normal vector 

defines the linear subspace 9 tangent to the collision hypersurface at q. 
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Fig. 6 schematically depicts the geometry of an elastic collision. If u(-) is the velocity before the collision, the 
velocity after the collision is such that the component in the plane 9 tangent to the hypersurface at q remains 
unchanged while the normal component changes sign, i.e., 

u (+) =u (-) _2( n . u (-)) n (112) 

which is the rule of elastic collisions (and the basic equation of geometrical optics) . 

For the following discussion we also need to define the linear subspaces 9 M which are perpendicular to the velocities 
uM. Following Sinai, we introduce the transformations: 

(a) The projection of Q^" 1 "' on to 3K~) , parallel to n, 

U : -» , given by U = 1 - 2 n n T , (113) 

which takes a vector in the hyperplane and projects it onto the plane such that the difference between the 
two vectors is parallel to n, and that the components of the two vectors in the directions perpendicular to n are equal. 
Here the superscript T denotes a transposed vector, or operator, and 1 is a unit operator. The transformation U is 
an invertible, orthogonal transformation given by an orthogonal matrix in the Nf dimensional configuration space 

U = U T = U 1 = 1 - 2 n n T , (114) 

with determinant det U = — 1. 

(b) The projection of onto 9 parallel to u'~) 

V : -> 9, 

is given by the matrix 

u(") n T 

v = 1 " • ^ 

(c) The projection of 9 onto S3( _ ) , parallel to n 

V T : 9 -> (116) 
is given by the transpose of the matrix V, Eq. (115). 

B. The Second Fundamental Form 

The convexity of the hypersurface <r(q) = determines the defocusing character of the collisions. The convexity is 
characterized by an operator known as the second fundamental form, 4 and which gives the variation of the normal 
vector with respect to a variation of a point on the hypersurface. That is, we define the second fundamental form 
K(q)by 

n(q+<kO - n(q) = K(qHq (117) 

where both q and q+^q belong to the hypersurface. Consequently, the variation Sq is perpendicular to the normal, 
n • Sq =0. Similarly, n 2 = 1 implies that n-<Sn = 0. As a result, the second fundamental form maps vectors, Sq, in the 
tangent plane onto vectors Sn in the tangent plane, 

K:3->3, (118) 

which is symmetric. If the operator is nonnegative, K > 0, the collisions are defocusing or neutral, and they may 
lead to a dynamical instability. In the case that K < the collisions are focusing or neutral, but, as illustrated by 



4 The first fundamental form determines the Riemannian distance between two nearby points on the surface. 
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Bunimovich's stadium billiard, the collisions do not necessarily lead to a dynamically stable situation. In order to see 
the meaning of the sign of the second fundamental form consider the inner product 

Sq T ■ Sn = (5q T • K • iq . (119) 

For surfaces that are defocusing, this inner product satisfies <5q T <5n >0, while for focusing surfaces <5q T <5n < as 
may be immediately checked by considering a sphere in three dimensions and constructing this inner product for both 
outward (Sq 1 ■ Sn > 0) and inward ( Sq 1 ■ Sn < ) normals. In the case of outward normals, 




where d is the radius of the sphere, lj_ = e#ej + e^ej is the unit operator in the tangent plane to the sphere at point 
q, and eg and are mutually orthogonal unit vectors in this tangent plane. 



C. An Example 

We now construct the various geometric quantities that describe a collision between particles 1 and 2. The hyper- 
surface of collision is locally defined by 

|n-r 2 | 2 = (n - x 2 f + ( Vl - y 2 f + ( Zl - z 2 f = d 2 . (120) 

A parametric representation of this hypersurface is given in terms of center of mass and relative coordinates by 

q= (Rl2 + ^€12, Rl2 - ^e 12 ,T 3 ,T 4 ,...,T N ) (121) 

where R = (r 1 +r 2 )/2 is the location of the center of mass of the two spheres 1 and 2 at the instant of the collision 
and £12 is the unit vector in the direction of the line joining their centers at the time of the collision. The unit vector 
which is normal to the hypersurface at the position Eq. (121) is given by 

n=-J=(£ 12 ,-£i2,0,...,0). (122) 

After applying the collision rule, Eq. (112) we find the usual result 

(123) 







-) _ 


[6l 2 -(v(->-v(->)]6 12 




= V2" 


-> + 


[ei2-(v(- ) -v(- ) )]6 12 




=v 3 - 


-) 





v (+) _ v (-) 

which conserves the total energy and momentum. 

The defocusing character of hard-sphere collisions is shown by considering the second fundamental form. The 
variation in position on the hypersurface is achieved by making a change in the unit vector £12 and in r 3 , rjy. Thus 
Sq is given by 

Sq= (^Se 12 ,-^Se 12 ,Sr 3 , ...,3r N ) , (124) 

and Sn is 

Sn=-j=(Se 12 ,-Se 12 ,0,...,0) . (125) 
From these two results it follows immediately that the second fundamental form is given as 

K = ^(1,1,0,0,...,0) . (126) 

Clearly the collision between the two spheres is defocusing as expected. However, defocusing only occurs in two 
directions while the other directions are neutral so that such billiards are referred to as semidispersing billiards [14]. 
One can think of the Nf dimensional configuration space for hard spheres as being bounded by hypercylinders with 
both defocusing and neutral directions much like ordinary three dimensional cylinders. 
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D. Linear Stability 



To develop a method to compute dynamical quantities for a hard-sphere gas we need to consider infinitesimal 
perturbations of some reference trajectory. The perturbations are best represented in a local frame of coordinates 
with one axis parallel to the velocity u 

Sq =6q\\e\\ + <5<?±ie ± i + Sq± 2 e ±2 H Sq± N f- 1 e ±Nf - 1 (127) 

where e|| = u with u • u =1. An equation similar to Eq. (127) obtains for the infinitesimal variation in velocity with 
respect to the reference trajectory. We remark that the condition u 2 = 1 implies that the velocity perturbation is 
necessarily perpendicular to u: u<5u = 0. Therefore we have 6u\\ = 0. Similarly we can set 6q\\ = 0, otherwise the 
perturbation can be assigned to another position along the trajectory. Consequently, the infinitesimal perturbation 
belongs to a linear space of dimension 2(Nf — 1). 

Our purpose is to obtain the time evolution of the infinitesimal perturbation over a time interval T which will define 
the monodromy matrix M(T) 

(&$) = ^ ■(£$)■ 

The time evolution is composed of collisions and free flights between the collisions. We shall first determine the 
monodromy matrix for a free flight and then the matrix for a collision. 



E. Free Flight 



During a free flight the velocity vector is constant although the position changes according to 

q(t)= u(0) t + q(0) , (129) 

where [q(0),u(0)] are the position and velocity after the previous collision. As a direct consequence, before the next 
collision 

where t is the time of flight between two consecutive collisions, which is also equal to the distance between impact 
points since the velocity u is normalized to unity. The monodromy matrix for a free flight is therefore 

Mf re e flight = ( q j j • (131) 



F. Collision 



If the collision illustrated in Fig. 6 is perturbed, the position and velocity just before the collision are displaced by 
(Jqj^\ Juj^) which both belong to the subspace perpendicular to the incident velocity u(~). As a consequence, 
the collision of the perturbed trajectory does not occur at the collision point q of the reference trajectory but at the 
nearby point q+^q on the hypersurface. we note that the perturbation Sq of the impact point belongs to the tangent 
subspace 9 and is determined from the perturbed position Sq^ by the projection from 9?( - ) onto 9 parallel to the 
incident velocity u(~', that is by Eq. ( 115) 

8q = V-dq^. (132) 

The perturbed trajectory after the collision is issued from the point q+^q with a velocity u( + )+5u( + ) given by the 
collison rule Eq. (112). The intersection of this trajectory with the tangent space 9 defines the perturbation in position 

Sq^ of the outgoing trajectory. Therefore the outgoing perturbation is given by the projection of the perturbation 
Eq. (132) of the impact point from the tangent subspace 9 onto the subspace parallel to the outgoing velocity 
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u'" 1 ") . In Fig. 6 we observe that the composition of the two successive projections, 9^ ) — > 9 — > is the projection 
from O^ - ' onto parallel to the normal vector n which is the inverse U 1 of the isometry Eqs. (113)-(114), 

<Jq<+> = U 1 -Sq^l (133) 

At the perturbed impact point, the normal vector is no longer identical to n but is perturbed as determined by the 
second fundamental form, Eq. ( 118). Accordingly, the velocity perturbation after the collision can be obtained in a 
straightforward way by differentiating the collision rule Eq. (112) 

5u { + ] = Su^ - 2(n-W L ~ ) )n - 2(n • u ( ~ } )K • 6q - 2(u^' • K • <5q)n (134) 

where <5q is the perturbation of the impact point, Eq. (132). Defining the angle <f> by 

cos <j> = n • u( + )= — n • u^ - ) (135) 

and using Eqs. (114) and (132), we find that Eq. (134) becomes 

5u ( + } = U" 1 -[<5uf L " ) +(2cos0) V T -K ■ V-Jq^] . (136) 
According to Eqs. (133) and (136), the monodromy matrix of a collision is 

M collision = ( (2cos0)u U -i. v t. k v J i ) (137) 
u 1 \ { 1 o 



U 1 J V (2co S(? !))V T -K • V 1 



G. Expanding and Contracting Horospheres 

The local Lyapunov exponents are determined by the rate of dispersion of trajectories in the vicinity of the reference 
trajectory. The dispersion is given by the "horosphere" which is a local sphere tangent to a front of trajectories 
accompanying the reference trajectory and issued from a common initial position in the past. The front is expanding 
so that we talk about the expanding horosphere which is nothing else than the local unstable manifold. For the case 
of one particle moving is a fixed array of scatterers this expanding horosphere is easy to describe, and it has been 
discussed in detail in Refs. [14,17]. 

The expanding front has a local curvature which is characterized by the second fundamental operator B u defined 

by 

<W = ^q± = B u -<Scll. (138) 
Accordingly, the local unstable manifold has the following parametric representation 

(&) = (bX) (139 » 

in the linear subspace. A similar representation holds for the local stable manifold. 

Let us consider a trajectory from an initial condition X = (q, u). Fig. 7 depicts the backward and forward portions 
of it. Collisions occur at the impact points n and times t n with n € Z. We denote by r n +i >n the time between the 
collisions n and n+ 1. Furthermore, B^~^(n) denotes the second fundamental operator of the horosphere immediately 
before the n th collision, i. e., at the end of the preceding free flight between the collisions n and n — 1. On the other 
hand, B^~\n) denotes the fundamental operator immediately after the n th collision. B u (t) denotes the fundamental 
operator during a free flight. 

We assume that the second fundamental operator of the expanding horosphere is fixed at some collision in the 
remote past and look for the operator at the initial condition. The operator is successively modified by the free flight 
and collisions according to the monodromy matrices Eqs. (131) and (137). The second fundamental operator is 
related to the one before the monodromy matrix M by 
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&) = (b&) = "-(bX) <-» 

which is solved by eliminating Sq± and 5q' ± between both lines. 

Applying this equation to the monodromy matrix Eq. (131) for a free flight , we find 

B^(n + 1) = [r n+1 ,„l + Bi+'fn)- 1 ] ~\ (free flight). (141) 

For a collision, we have from Eq. (137) 

B( 1 + )(n) = U; 1 -[(2cos0„)V^K n -V„ + B( 1 -)(n)]-U„ (collision). (142) 

Combining Eqs. (141) and (142) for successive backward collisions, we obtain Sinai's matrix continued fraction 
expression for B u (t): 



B u (t) 



rl + U- 1 • [(2cos&,X-K„.V n + = ]-U„ 



T„,„_il + U„i! • ^COS^n-OV^.j-Kn-i-Vn-i + — ] U„_i 

Tn-l,n-2* + " " " 

(143) 

with t = t — t„ and for t n+ i > t > t„. Because K n > 0, and cos(j> n > 0, the matrix B u (t) > so that the expanding 
character is maintained during the whole time evolution in the case of a hard-sphere gas. Let us emphasize that 
the operator (143) is defined locally for each initial condition X and can be obtained by integrating backward the 
trajectory $'X for — oo < t < to determine the successive past collisions and the corresponding quantities appearing 
in Eq. (143). 

To see the connection between the second fundamental operator and the Lyapunov exponents of the hard-sphere 
gas we proceed as follows. We use the fact that between collisions the quantity Sq± develops as 

8<l±(t n ) = <5qj_(£„_i +r n> „_i) = 5q±(t n - 1 ) + T„ in _i<Suj_(* n _i) = [1 + r„ in _iB^ +) (n - 1)] -<5q±(i„-i) (144) 

where the quantities are the values immediately after the collisions, as given by Eq. (142). After a sequence of n 
collisions, Sq±(t) is given by 

<fqx(t n )= [l + T„, n _ 1 Bi + '(n-l)]-[l + T„_ 1 ,„_ 2 B( I + )(n-2)]-.-[l + T 1 , Bi + '(0)]- ( 5q ± (0) , (145) 

with to = 0. If there is an exponential separation of trajectories we would expect that for long times ||<Sqj_(T)|| ss 
(exp AT) 1 1 5q(0) 1 1 . The exponential growth factors A would satisfy the relation 



VA,= lim ^lndet[l + r„,„_ 1 B( I + '(n-l)]---[l+T lo B( I + )(0)] (146) 

Ai>0 

1 n— 1 

lim -^lndet[l + r i+M Bl+)W] 

8=0 

n— 1 

lim i / T,+1 " drtr [rl + BW(i)- 1 ] -1 (147) 



T^oo T 

Now by using Eq. (141) in the form 

B u (t n + r) = B(+)(n)[l + rB( + )(n)] 

we readily find that 



1 f T 

A< = lim - / dttrB u (t) . (148) 
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In the event that the system is closed and ergodic, then the time average in Eq. (148) can be replaced by an ensemble 
average with respect to the equilibrium invariant measure of the system, so that in this case 



^2 Aj = // eq (tr B„) 



(149) 



A;>0 



for almost all trajectories. For open systems, a similar result obtains but the appropriate nonequilibrium measure 
must be used in computing the average of the trace of the second fundamental operator. We can therefore identify 
the local dispersion rate for billiards as 



The form Eq. (150) plays for hard spheres the role of Eq. (32) for systems with smooth potentials. With 
respect to smooth Hamiltonian systems, a simplification for billiard systems comes from the fact that the local 
Lyapunov exponents are given directly in terms of quantities which can be constructed from successive collisions. It is 
possible that numerical calculations of the sum of positive Lyapunov exponents for billiard systems could be calculated 
efficiently and quickly using Eqs. (146) or (148). It is worth noting that the individual Lyapunov exponents may also 
be calculated by using the second fundamental operator (143). 

For the special case of a billiard system consisting of one moving particle in a regular triangular array of fixed hard 
disk scatterers - the triangular Lorentz gas - one can evaluate the continued fraction, Eq. (143), in terms of scalar 
quantities, and the positive Lyapunov exponent has been determined numerically for both open and closed systems 
by Gaspard and Baras [16,17]. For open systems, the positive Lyapunov exponent describes the expanding manifold 
for the trajectories on the fractal repeller. In the case of random Lorentz gases - where a particle moves in a fixed but 
random array of hard disk or hard-sphere scatterers - van Beijeren, Dorfman, and Latz [18,44] have determined the 
sum of the positive Lyapunov exponents for both closed and open systems in the case that the density of scatterers 
is low. 

As a consequence of the considerations presented here we see that even for billiard systems it is possible to define the 
sum of the positive Lyapunov exponents in terms of the average value of the trace of the second fundamental operator. 
For closed systems, this quantity determines the Kolmogorov-Sinai entropy of the billiard, by Pesin's theorem. For 
open systems, this quantity provides the sum of the positive Lyapunov exponents for trajectories on the repeller, 
which is an essential ingredient needed for the escape-rate formalism for open systems. 



As we remarked in subsection II. F and as is illustrated in Fig. 5, the hard-sphere gas with periodic boundary 
conditions or in a rectangular box presents special periodic orbits for which the particles do not undergo collisions. 
As a result the second fundamental operator Eq. (143) is of free-flight type, Eq. (141) which decays to zero as t~ x as 
t ^ oo because K = for the particular case of collisions on flat walls. Therefore, all of the Lyapunov exponents of 
those special orbits vanish and the system is nonhyperbolic. These special periodic orbits form continuous families of 
dimension Nf — 1. We discussed in subsection II. E how the large-deviation formalism can be applied to such systems. 

Since the transport coefficients are bulk properties obtained after a thermodynamic limit, the theory of paper I is 
also valid for a system which is modified at its boundaries. Using the independence of the transport coefficients of 
surface effects, we can consider a hard-sphere gas in a box with convex walls as in Fig. 8. Fig. 8a shows walls given 
by portions of spheres with a radius of order of the box size. In this case, the convexity is global and mild but enough 
to turn the system into a hyperbolic system. Indeed, the families of special orbits have now disappeared and all of 
the periodic orbits are now unstable. 

Another possibility is shown in Fig. 8b where the walls are composed of many portions of small spheres modeling 
the atoms of the walls. Here also, the continuous families of special periodic orbits have disappeared. There may still 
remain rare periodic orbits where the particles have no mutual collisions but these periodic orbits are now unstable 
since the walls are defocusing. Therefore, the system here is also hyperbolic and the escape rate formula Eq. (100) 
strictly applies. 




(150) 



A;>0 



H. Marginally Unstable Periodic Orbits 
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IV. LATTICE-GAS AUTOMATA 



A. Lattice-gas automata as Markov chains 

Lattice-gas automata are probabilistic Markov chains over system states defined by the set of discrete states taken 
by individual particles. The single particle state for particle takes one among m integers which determine its 
position and velocity. The collection of all the single-particle states (£i,£2, £3, • • • , £zv) = w defines a state of the whole 
system. As a consequence the total number of possible system states is m N [45]. This set of possible system states 
plays the role of an alphabet A = {l, 2, 3, . . . , m N } of a symbolic dynamics M = £.4 in which the trajectories of the 
system correspond to bi-infinite sequences of states 

W = . . . . . . . (151) 

The successive symbols u>k give the state of the system at successive times ty. = kAt. The set of all bi-infinite sequences 
defines the phase space of the lattice-gas automaton. This so-defined phase space is a continuum, as it should be. 
For example a simple automaton with an alphabet of two symbols, and 1, and a set of trajectories consisting of all 
bi-infinite sequences of these two sysmbols can be mapped onto a unit square, as is done in the baker's transformation. 
We remark that a symbolic dynamics establishes a one-to-one correspondence between the system trajectories and 
the bi-infinite sequences. 

Gaspard and Wang have shown elsewhere [46] that an area-preserving map can be constructed which is isomorphic to 
the Markov chain. In this way we can establish a connection between some properties of the Markov chain and those of 
a deterministic dynamical system. However, the dimensionality is lost in this connection. As a consequence, geometric 
properties of dimension which are typical of differentiable dynamical systems cannot be recovered. Nevertheless, we 
are here interested in the quantities appearing in the escape-rate formalism. In particular, a formula like Eq. (2) 
only contains the sum of positive Lyapunov exponents, a global quantity which does not require the knowledge of 
individual Lyapunov exponents. Therefore the escape-rate formalism can be extended to lattice-gas automata. 

The lattice- gas automaton is fully characterized by the matrix P of the Markov chain on the system states. The 
elements P uu i of the matrix P give the probabilities of transition between two successive states to and to'. Since 
probability is conserved and the Markov chain admits an invariant vector we know that 

^-Puu'=l, and ^PuPuu' = Pu' (152) 

with J2 u Pu> = 1- We may assume that the Markov chain is ergodic so that the invariant vector is unique. The 
invariant probability measure is then defined as 

(153) 



B. Escape- rate formalism of Markov chains 

Although the dynamics takes place within the full set of states, A, there may exist a subset of states, B C A , 
which are visited in a transient way. This might happen, for instance, if we impose a nonequilibrium constraint on the 
system, such as requiring that a Helfand moment lie in an interval Eq. (1). If the initial state belongs to this subset, 
the state will escape from this subset after a finite time for almost all of the trajectories Eq. (151). Nevertheless there 
exist trajectories Eq. (151) which remain forever on this subset B. These trapped trajectories are exceptional in the 
sense that their measure, Eq. (153), is vanishing and consequently they form a repeller TZb = Eg- The first time 
at which a given trajectory escapes from the subset B defines a problem of first passage. The dynamics of the full 
system with respect to the repeller is similar to a scattering process where the trajectory visits for a while a vicinity 
of the repeller. In this image where a trajectory is going in and out of the subset B , the repeller may be considered 
as a predefined nonequilibrium fluctuation, as discussed elsewhere [47,48]. 

Since the dynamics is transient on the repeller there exists an escape rate 7. We can calculate this escape rate by 
constructing the Markov subchain between states of the subset B. The values of the transition probabilities are given 
by a submatrix which is contained in the matrix P of the full Markov chain Eq. (152). This submatrix Q is composed 
of the elements of P between the states of the subset B 

A\B B 

_ A\B ( * * \ _ 

P = B { * Q J ' < 154 ) 
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where * denotes other submatrices. 

Because it is obtained by truncating the full matrix P, the submatrix Q is not a stochastic matrix obeying Eq. 
(152). In particular, the leading eigenvalue of Q is no longer 1 but is smaller than 1. Actually the eigenvalues of the 
submatrix Q give the decay rates out of the subset B so that the leading eigenvalue gives the escape rate 7 according 
to 

Q\v) = exp(— "f)\v), and (u\Q = exp(— ^)(u\ (155) 

which defines the leading right and left eigenvectors of Q. The escape rate is here in units of the time steps, At, of 
the automaton. If Q is a nonnegative, irreducible, aperiodic matrix then according to the Perron-Frobenius theorem 
all of the other eigenvalues are less in absolute value the the largest eigenvalue, exp(— 7). 

Using the right and left eigenvectors we can construct the Markov subchain of the repeller dynamics in terms of 
the following matrix and its associated invariant vector 

n wJ = exp(7)Q^'— , and tt u = (156) 

Vui {u\ v) 

which is now stochastic since 

^ /Lw' =1, 7r "-^"' = n u'> an d ^ 7r M = 1. (157) 

oj'EB gjEB gjEB 

This stochastic matrix defines an invariant measure on the repeller which is given by 

fi B (uio^l ■ ■ ■ w„_i) 

— ^uo-^uouin ulU2 . . . n UJri _ 2UJri _ 1 (158) 

where u!k G B. This new invariant measure /xg gives nonvanishing probabilities only for trajectories staying on the 
repeller. 



C. Lattice-gas automata as chaotic dynamical systems 

We now have the necessary elements to proceed with a derivation of the escape-rate formula. The escape rate 7 
was already defined in terms of the leading eigenvalue of the submatrix Q. On the other hand, the KS entropy of the 
dynamics on the repeller is defined as the KS entropy per time step At of the Markov chain Eqs. (157) and (158), 

/iks(/«b) = - ^2 v un uu >\nII uul . (159) 

We need to find the quantity which plays the role of the sum of the positive Lyapunov exponents, X^a >o ^ ls 
not obvious at first sight how to do this. The solution can be found by mapping a Markov chain onto a deterministic 
map [46,49], and then computing the sum of the positive Lyapunov exponents for the deterministic map. One finds 
in this way that the inverses of the probabilities, Q U u' , for the separate steps of the Markov chain play the role 
of the stretching factors by which trajectories separate in the map or by which the probability is dispersed in the 
Markov chain. The inverse of the dispersion factor for the trajectories visiting successively the states loqU}\loi . . . Lo n -i 
is therefore given by 

eXp-[/(w Wl . . .U) n -l) = Qu UlQuiU2 ■ ■■Qu n -2U n - 1 - (160) 

The role of the sum of the Lyapunov exponents is then played by the average quantity 

Hb{u)= lim - S~] /xb(w . . . w„-i) lnexp[/(w . . .cj„_i), (161) 

n— >oo fi — * 

uo...u„_i£B 

with the definition Eq. (160). Eqs. (160)-(161) define for lattice-gas automata the analogue of the dispersion rates 
given by Eq. (32) for systems with smooth potentials and by Eq. (150) for hard spheres. 

Using the factorization property of the dispersion factors, Eq. (160), and the Markov property of the measure Eqs. 
(157) and (158), we obtain 

Hb{u) = ^2 K^n^u' In 1 (162) 
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which is positive since the matrix elements are smaller than one: Q^^i < 1. The dispersion rate u has the same 
units, Ai _1 , as the escape rate 7 and the KS entropy, Eq. (159). The isomorphism between the Markov chain and 
a deterministic, area-preserving map has been discussed in detail for the multibaker transformation and for Lorentz 
lattice gases [48,49]. This isomorphism shows that the sum of the positive Lyapunov exponents for the deterministic 
systems is identical to /xs(u) given above, 



Hb{u) 



E 

Ai>0 



(163) 



both calculated on the repeller Eg. This last property guarantees that the averaged quantity u is the unique analog 
of the sum of the positive mean Lyapunov exponents for lattice-gas automata. 

Replacing the definitions Eq. (156) of the matrix II and of the invariant vector 7r in the expression Eq. (159) for 
the KS entropy, we obtain the identity 



7 = hb(u) - h K s(pB) 



(164) 



which is the escape-rate formula (100) for lattice-gas automata. This formula gives the escape rate from the repeller 
Eg as the difference between an average dispersion rate playing the role of the sum of the positive Lyapunov exponents 
and the KS entropy. 

Similarly, we have the analog of the Pesin theorem. When we relax the constraint on the subset of allowed states, 
then A = B, the repeller becomes the full phase space, and Q=P. As a consequence the escape rate vanishes and 



(i A (u) = h KS (/j, A ). 

We can also define the dynamical pressure function of section II for lattice-gas automata according to 



P(d) 



lim — In 

n— ¥00 fl 



E 







(165) 



(166) 



This pressure function has all of the properties described in section II, which relate it to the preceding quantities like 
the escape rate, the dispersion rate u and the KS entropy 



7 



dP 



and /iKs(jUs) 



/3=i 



(167) 



0=1 



It is a straightforward calculation to show that the quantities 7, hb(u), and /iks(A'b) derived from the pressure 
function are identical with those given by Eqs. (155), (162), and (159), respectively. This is accomplished by using 
the spectral decomposition of the matrix Q, in terms of its eigenvalues and eigenvectors. If this matrix is irreducible 
and aperiodic the dominant behavior of Q", for large n is determined the the largest eigenvalue and the associated 
left and right eigenvectors. The correspondence of the two kinds of expressions for the dynamical quantities then 
follows immediately. 

Moreover, the analog of the topological entropy per time step At is given by 



^top(B) = P(0). 



(168) 



The topological entropy is the largest eigenvalue of the topological transition matrix associated with the Markov 
subchain II, i.e., the matrix with elements or 1 depending upon whether n uu < is zero or positive. If the subchain 
matrix Q is a strictly positive Mg x Mg matrix with M& < m N , then the topological entropy is equal to h top (B) = 
In M B < iVlnm. 



V. RELATIONSHIPS TO OTHER THEORIES AND APPROACHES 



Recently, there has been substantial progress in understanding the relationship and connections between a number 
of approaches for describing transport processes in classical systems on the basis of dynamical systems theory. In this 
section, we briefly describe this progress and indicate some directions for further work. 
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A. Spectral theory of the Perron- Frobenius operator 



In classical statistical mechanics, the time evolution of the probability density p is described by the Liouville 
equation 

d t p = L p where L = {H, ■ } (169) 

is the Poisson bracket with the Hamiltonian. The solution of Eq. (169) is formally given by pt = exp(Lt)po. 
Alternatively, the time evolution on an energy shell Ai is described by the Perron-Frobenius operator 

MX) = M^'X) = / «S(X-$'Y) po(Y) dY =(P'po)(X) (170) 

J M 

This operator has a kernel which is a Schwartz distribution rather than a function, indicative, of course, of the 
unique time development of a given phase point, Y. This is to be contrasted with the type of kernels appearing in 
the expressions for the time development of probability densities that characterize stochastic processes, such as in 
solutions to the Fokker-Planck equation, which involve a standard Green's function. Resolving this difference, that is, 
showing how deterministic systems can show stochastic-like behavior, is at the heart of the major difficult problems 
of nonequilibrium statistical mechanics. Recently, some of these difficulties seem to have been overcome with the 
development of a spectral theory of Perron-Frobenius operators for hyperbolic systems, in particular, by Pollicott 
and Ruelle [50]. These authors have shown that the Perron-Frobenius operator admits resonances corresponding to 
exponentially decaying eigenstates 

P l \K) = e s "' \4>n) for t >0 (171) 

The eigenstates {|0n)} are distributions without any corresponding density which explains the failure of previous 
attempts to define decaying eigenmodes in terms of density functions. The resonances characterize the forward semi- 
group so that Res„ < 0. By time reversibility, antiresonances are associated to each one of them, which characterize 
the backward semigroup. The forward time evolution of an observable A can be decomposed like 

(A|p'|po) = E < A \M eSnt $M + ■■■ ( 172 ) 

n 

where the dots denote other terms due to a possible Jordan-block structure [51-53]. The eigenstates {|0n)} and their 
adjoint {{0n|} are distributions which are defined on sufficiently smooth observables A and initial densities po- 

For closed systems which are ergodic and mixing, the leading resonance is vanishing Sq = and it corresponds to 
the unique microcanonical invariant measure so that (A\(f> ) = n e q{A) and (4>o\po) = A'eq(po)- The higher resonances 
give the relaxation rates toward equilibrium in the forward semigroup. 

For scattering systems which are decaying, the leading resonance gives the escape rate Sq = —7 [54,53]. The 
corresponding eigenstate \<f>a) is concentrated on the unstable manifolds of the repeller: Cl[W / tt (7?.)]. Such a state 
corresponds to the construction of the fractal by the process Eq. (80) and the definition of the escape rate by Eq. 
(88). By time reversibility, the values of the escape rate obtained in both theories are identical according to Eqs. 
(87)-(89). 

In the context of transport phenomena, one considers a spatially extended system, which often can be suitably 
described by imposing appropriate periodic boundary conditions [53,55]. The dynamics over the full phase space can 
then be reduced to the dynamics inside a basic cell of the system by spatial Fourier transforms which introduce the 
wavenumber k. The density can be decomposed into components pk of wavenumber k which are quasiperiodic in 
space according to Bloch's theorem. These components evolve in time under the following Perron-Frobenius operator 
which now depends on the wavenumber k [56] 

(P*p k )(X) = J e ik a < x ) S(X - $'Y) p k (Y) dY (173) 

where a(X) is a function from the phase space to the periodic /-dimensional physical space which describes the jumps 
of the trajectories at the periodic boundaries. This Perron-Frobenius operator also admits Pollicott-Ruelle resonances 
and corresponding eigenstates 

Pt\4>kn) = e s " (k)i \fc n ) for t > (174) 
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depending on k. In the case of diffusion, the leading resonance so(k) gives the dispersion relation: so(k) = —Dk 2 + 
0(k i ) [53,56]. This result shows the connection with the chaotic scattering approach were the leading Pollicott- 
Ruelle resonance or escape rate is so = — 7 = —D(tt/L) 2 + 0(L~ 3 ) so that we can identify the wavenumber of the 
corresponding hydrodynamical mode as k = (tt/L) + (D(L~ 2 ) [2]. 

For the multibaker model, the connection between the Pollicott-Ruelle resonances of the closed and scattering 
systems have been described in detail in Ref. [48] and the spectral decomposition Eq. (174) has been explicitly 
constructed in Ref. [53]. Tasaki and Gaspard [57] have further shown that a nonequilibrium steady state corresponding 
to a constant concentration gradient is obtained in the zero- wavenumber limit as c^kl^ko) | k _ of the leading eigenstate 
of Eq. (174). Accordingly, such nonequilibrium steady states are also given by distributions. Moreover, Tasaki 
and Gaspard have proved for the multibaker model that the nonequilibrium steady states are described by singular 
measures having for support the full phase space and, nevertheless, exhibiting self-similar properties. A remarkable 
fact is that these self-similar properties have their origin in the fractal properties of the repeller and, in particular, 
of the escape-time functions of Eqs. (76) and (77), as shown by taking the large-system limit (B — > Ai) in order to 
connect the open and scattering configuration to the closed and infinitely extended configuration. 



B. Periodic-orbit theory 

By using periodic-orbit theory [56,58,59], one can obtain in hyperbolic systems the characteristic equation for the 
resonances of the Perron-Frobenius operator in terms of the periodic orbits which thus appear to be directly related 
to the kinetics of the system. The periodic-orbit theory is based on the trace of the Perron-Frobenius operator which 
is given by 

p r—1 ^ 

where the first sum extends over the prime periodic orbits {p} and the second sum over their repetitions. T p is the 
prime period of p, M p the linearized Poincare mapping in a surface of section transverse to the periodic orbit, a p is a 
spatial vector representing the distance travelled in the full phase space when the orbit closes on itself in the reduced 
phase space. The Laplace transform of this trace can be expressed as 



I 



00 - 8 

e- st tr P£dt = — In Z(s) (176) 



ds 

in terms of a Selberg-Smale zeta function defined by a product over the prime periodic orbits 

(m 1 +l)...(m i +l) 



*w = n n 

p m i . . . m, l = 



_ exp(zk-a p - sT p ) 

nLiAi°iAi° m< 



(177) 



Here the quantities A p ' are the expanding eigenvalues of the matrix M p , and L is the total number of expanding 
directions. The zeros of the Selberg-Smale zeta function are the Pollicott-Ruelle resonances; Z[s n (k)] — 0. Zeros are 
obtained after transforming the infinite product into a series by cycle expansion. In this way, a formula has been 
obtained [56] which expresses the diffusion coefficient in terms of the Lyapunov exponents of each periodic orbit, their 
periods, and the vectors a p 

n = 1 EgLo E P1 ... P „ (-)"(a pi + • ; • + a p „ ) 2 /|A pl • . • V| 

4 E~ o E P1 ... P „ (~) n (T Pl + ■ ■ ■ + T Pn ) V|A P1 • • • A p J [ > 

for a two-degree-of- freedom system for which there is a single unstable direction (L = 1). 

The periodic-orbit method has been applied not only to the study of diffusion in the Lorentz gas [56] and in other 
models [60] , but also to the calculation of the Pollicott-Ruelle resonances of the three-disk scatterers [54] as well as 
of open and closed configurations of the multibaker [48] . The method is particularly interesting to obtain resonances 
beyond the escape rate which are difficult to obtain by direct simulation methods. However, the method is very 
demanding because a large number of periodic orbits may be required as it is the case in the Lorentz gas. 

We conclude that there exist direct connections between the theory described in the present paper and the periodic 
orbit theory of the Perron-Frobenius operator. In both theories, the very same escape rate can be obtained as well as 
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the other characteristic quantities of chaos such as the Lyapunov exponents and the partial Hausdorff dimension, as 
described elsewhere [16,17]. 

Recently Morriss, Rondoni, and Cohen have applied periodic-orbit theory to obtain the equilibrium pressure exerted 
by moving particles on the scatterers in a periodic, triangular Lorentz gas [61]. By carrying out computer simulations 
these authors were able to show that within the numerical errors, the equilibrium-ensemble measures are equivalent 
to dynamical measures based upon the stretching factors Ap ' whose logarithms are the positive Lyapunov exponents 
for the corresponding periodic orbit. This result provides a further support to the ideas of Sinai, Ruelle, and Bowen 
on the dynamical foundations of the equilibrium measures. 

C. Thermostatted-system approach 

In this approach, the system is subjected to some external field such as an electric field and to a special ther- 
mostatting force which maintains a constant kinetic energy [62-67]. These systems are deterministic and the dy- 
namics is given by a set of ordinary differential equations, X = F(X), which is time-reversal symmetric under 0: 
(q, p, t) — > (q, — p, —t). However, contrary to Hamiltonian systems where phase-space volumes are preserved divF = 0, 
the dynamics of thermostatted systems do not preserve phase-space volume but 

div F < for X e M+ 

div F > for X e M- (179) 

so that the phase space is composed of two parts Ai = Ai+ U M- which are mapped onto each other by time reversal 
OAi+ = Ai-. The divergence of the vector field appears to be directly related to the transport properties by the 
way the thermostatting force is constructed. Moreover, because of the assumed coupling to a thermostat, an entropy 
production is defined in this approach. 

As a consequence of Eq. (179), the domain Ai+ attracts the future of the trajectories and contains an attractor A+. 
By time reversibility, the other domain Ai- contains a repeller A- which is repelling in every directions (contrary 
to the repellers 1Z considered in the present paper which are of saddle type). Under certain conditions, the attractor 
A+ may be a chaotic and fractal attractor to which the thermodynamic formalism described here is of application 
as shown by Chernov et al. [66]. These authors show, among othert things, that the attractor is the support of a 
Sinai-Bowen-Ruelle invariant measure for which the Pesin formula holds, and the escape rate vanishes, 7 = 0. 

Recently, Gallavotti and Cohen [68] have described in some detail the construction of the SRB measures for the 
attractor, using methods based on Markov partitions. In this approach, they were able to show that the SRB 
measures can be used to provide the theoretical foundations for, and to quantitively describe the numerically observed 
fluctuations of the thermostatted-system entropy production in a shearing fluid far from equilibrium. Here, too, one 
has a clear example of the utility of constructing dynamical measures to describe nonequilibrium phenomena. 

For thermostatted systems, the transport coefficients are related to the sum of all the Lyapunov exponents, A,, 
which is negative in the domain Ai+ and, therefore, on the attractor A+ because of Eq. (179), in contrast to 
our escape-rate formula Eq. (2) in which the transport coefficient is related to the difference between the sum of 
positive Lyapunov exponents and the KS entropy for trajectories on the repeller. Nevertheless, there appear to be 
important structural similarities between the two expressions for transport coefficients, which have been commented 
upon in the literature [69,70] since the two kinds of expressions are given by the ratios of dynamical quantities to 
the squares of external parameters - the system size in the escape-rate formula, and the external field strength in the 
thermostat formalism. In both approaches, the hydrodynamical scale is related to the kinetic and chaotic time scale by 
a mechanism of dimensional loss in phase space. However, the mechanism of dimensional loss turns out to be different. 
In the chaotic-scattering approach, the dimensional loss is due to the finiteness of the scatterer and to the escape 
resulting from the nonequilibrium boundary condition. On the other hand, the non- volume-preserving thermostatting 
force is at the origin of dimensional loss in the thermostatted-system approach. We might mention here also that a 
somewhat different approach to dimensional reduction in systems with shearing flows has been developed by Chernov 
and Lebowitz [71] based upon the construction of special boundary conditions which maintain the shearing flow. The 
connections between all of these methods still remain to be established. 

In section IV, we have seen that fractal dimensions can also be defined for probabilistic systems. Similarly, the 
concept of KS entropy has been generalized into a concept of e-entropy to characterize stochastic processes [46]. This 
observation shows that the concepts of fractal and of chaos are not restricted to deterministic systems. Accordingly, 
we can imagine that it may be possible to characterize a mechanism of dimensional loss even in stochastic systems 
which describe systems coupled to thermal baths and in which the bath degrees of freedom are taken into account 
by fluctuating forces along with dissipative forces as in Langevin processes. In this larger perspective, we think that 
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connections may exist between the descriptions of the transport and dynamical properties of thermostatted systems, 
and those with stochastic boundary conditions, with the chaotic-scattering approach described in this paper. 



VI. CONCLUSION 

In this paper, we have reviewed the derivation of the escape-rate formula for a variety of systems of interest for 
studies of transport processes in fluids. The derivation is based on the local stretching rates which we introduced in Eq. 
(30) for smooth hyperbolic systems. Beside smooth hyperbolic systems, we have also considered the billiard systems 
and the cellular-automata lattice gases. The relation between the escape-rate formalism and transport coefficients 
was discussed in I, and a number of examples of the use of this formalism have appeared in the literature recently 
[16-19,48]. It is of course essential to explore a number of specific cases in order to develop some intuition about the 
properties of the fractal repellers that control the escape-rate, and to study cases where transport is anomalous in 
some way in order to determine the limitations of this approach. 

There are other approaches to expressing transport coefficients in terms of the dynamical properties of the system 
under consideration. Such approaches have involved cycle expansions of dynamical zeta functions, or other properties 
of the zeta functions, the effects of Gaussian thermostats on the Lyapunov exponents of the system, as well as studies 
of the nature of the hydrodynamic modes that describe transport in the system. The connections between these 
various methods are only partially understood today. Clarifying these connections as well as establishing connections 
with more traditional approaches to transport theory based upon kinetic equations or Green-Kubo formulae remains 
a fruitful area for further study. 

For a deeper understanding of the escape-rate formalism, it is essential to find methods of determining the KS 
entropy for trajectories on the repeller in a way that does not rely on the escape-rate formula. For certain uniformly 
hyperbolic systems, this can be accomplished by taking advantage of the equality of the KS and topological entropies, 
since the topological entropy can be determined in terms of the eigenvalues of the topological transition matrix [48,72]. 
However, to have more general methods, we will need to develop a deeper understanding of the dynamics on the fractal 
repeller with a possible development of a theory of Markov partitions for such trajectories. Another promising approach 
is based on the "Maryland" algorithm [73]. Here, for two-degree-of-freedom systems, one calculates the Hausdorff 
codimension of the repeller by studying the intervals of continuity of the escape-time function such as that illustrated 
in Fig. 4. One can then show that the Hausdorff codimension is related to the KS entropy in the large-system limit. 
This method has been applied to the determination of the diffusion coefficient for the Lorentz gas on a triangular 
lattice [16,17]. However, the generalization of such considerations to higher-dimensional fractals is still in its infancy 
[74,75]. 

Finally, we mention again that the Sinai method for determining the dynamical properties of billiard systems can 
be used as the basis of a more rigorous kinetic theory of gases, based not on the assumption of molecular chaos, but 
based rather on the dynamical chaos of billiard systems. In this direction, Bunimovich and Spohn have been able to 
provide a rigorous proof of the existence and positivity for the shear and bulk viscosities of a periodic billiard fluid 
with two hard disks per unit cell [76]. Moreover, Simanyi and Szasz [77] have shown that it may indeed be possible 
to complete Sinai's arguments for a rigorous proof of the ergodic behavior for a system of hard spheres. 
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Figure captions 

Fig. 1. Schematic representation of the stable W S (K) and unstable W„(X) manifolds of the trajectory from initial 
condition X. £ is a surface of section transverse to the orbit. The invariant manifolds are shown on the respective 
sides of £ which are used for their construction. 

Fig. 2. Schematic behavior of the multivariate pressure function P(/3i,/3 2 ) for a three-degree-of-freedom system. 

Fig. 3. Schematic behavior of the pressure function in the cases of: (a) a closed hyperbolic system; (b) an open 
hyperbolic system; (c) a closed nonhyperbolic system; (d) an open nonhyperbolic system. 7 denotes the escape rate 
and 7eff are the effective escape rate; u the sum of the mean positive Lyapunov exponents; /iks the KS entropy per 
unit time; h top the topological entropy per unit time. 

Fig. 4. Typical behavior of the escape-time function versus the initial condition Xq, here for the ID logistic map 
x n +i = [ix n (l — x n ) with [i = 4.01. 

Fig. 5. Examples of marginally unstable periodic orbits with zero Lyapunov exponents in the cases of (a) the stadium 
billiard, (b) the Sinai billiard, (c) the hard-sphere gas with periodic boundary conditions, (d) the hard-sphere gas in 
a rectangular box. 

Fig. 6. Geometry of an elastic collision, n is the vector normal to the hypersurface of the billiard, are respectively 
the velocities before and after the collision. 9 is the tangent subspace. Sf^' are the subspace perpendicular to the 
ingoing and outgoing velocities. Sq^ are the infinitesimal perturbations in position before and after collision. Sq is 
the infinitesimal perturbation of the impact point q. The curvature of the collision hypersurface is not represented 
for the following reason. The hypersurface curvature plays a role only in the infinitesimal perturbations in velocities 
Su^ which are straightforwardly derived from the collision rule (112) by analysis as done in Eqs. (134)-(136). This 
analytical derivation does not particularly require a geometrical visualization. 

Fig. 7. Expanding and contracting horospheres of a billiard. These horospheres are geometrical representations of 
the local unstable and stable manifolds, wi' oc ^(X) and W / i' oc ^(X), at the phase-space point X = (q, u). (a) Backward 
trajectory determining the expanding horosphere. (b) Forward trajectory determining the contracting horosphere. 

Fig. 8. Hard-sphere gas in a box with convex, defocusing walls: (a) Each wall is composed of a portion of a single 
large sphere; (b) of the portions of many spheres modeling the atoms of the wall. 
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